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ABSTRACT 


The  NPS  Aerolab®  Low  Speed  Wind  Tunnel  located  in  Halligan  Hall  of  the 
Navy  Postgraduate  school  has  been  in  operation  since  1953.  Although  the  tunnel  is 
well  maintained,  its  data  acquisition  system  has  not  kept  pace  with  modem 
technology.  An  effective  but  affordable  solution  for  acquiring  data  was  needed.  It  was 
determined  that  a  software  package  known  as  Lab  VIEW©  would  provide  a  low  cost,  data 
acquisition  solution  that  will  enhance  the  capabilities  of  the  wind  tunnel,  while  at  the 
same  time  making  it  more  user  friendly  to  faculty  and  students.  The  focus  of  this  thesis  is 
the  design  of  a  VI  that  will  collect  and  plot  force  and  moment  data  from  a  six  component 
strain  gauge  balance  and  yield  real  time,  non-dimensional,  force  and  moment  coefficients 
in  six  degrees  of  freedom.  Wind  tunnel  tests  consisting  of  angle-of-attack  sweeps  in  the 
NPS  Aerolab®  low-speed  wind  tunnel  were  conducted  to  verify  (L/D)Isp  optimized,  M„ 
=  6,  conical-flow  waverider  data  obtained  in  1994  using  a  different  data  acquisition 
system.  Results  of  current  testing  substantiate  the  Lab  VIEW©  code  and  the  validity  of 
the  1994  test  data.  Analysis  of  the  current  wind  tunnel  test  data  resolved  pitching 
moment  concerns  related  to  the  1994  data. 
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I.  INTRODUCTION 


A.  THE  NEED  FOR  AN  IMPROVED  DATA  ACQUISITION  SYSTEM 

Perhaps  the  most  challenging  part  of  wind  tunnel  experimentation  is  acquiring 
data  that  can  be  easily  processed  and  analyzed.  The  wind  tunnel  data  acquisition  system 
in  use  prior  to  this  thesis  utilized  an  IBM©  PS2  computer  with  an  installed  data 
acquisition  card  for  obtaining  the  data.  In  this  manner  it  was  similar  to  a  typical  PC 
based  system,  but  in  other  areas  its  differences  were  readily  apparent. 

The  operating  system  of  the  IBM  PS2  is  not  compatible  with  today’s  PC, 
complicating  data  retrieval  into  existing  spreadsheet  and  computational  programs.  This 
incompatibility  carried  over  to  the  data  acquisition  program,  Microsoft  Quick  Basic,  used 
by  the  PS2.  Quick  Basic  is  a  programming  language  not  familiar  .to  most  faculty  and 
students  using  the  wind  tunnel  nor  typical  of  a  modem  data  acquisition  system.  The 
output  of  the  program  yielded  only  dimensional  forces  and  moments.  The  user  would 
then  have  to  non-dimensionalize  the  results  off-line  using  another  software  application, 
which  also  meant  that  real  time  results  could  not  be  displayed  during  the  data  acquisition 
process  in  order  to  compare  test  values  with  expected  results. 

Most  important  however,  the  PS2's  hard  drive  was  no  longer  functioning  (it  was 
broken).  With  the  above  factors  in  mind,  it  was  thought  more  beneficial  to  develop 
another  data  acquisition  system  for  the  wind  tunnel  rather  than  try  to  repair  the  existing 
one,  if  costs  could  be  kept  reasonable.  Developing  a  new  system  would  create  an 


l 


opportunity  to  update  data  acquisition  methods  for  the  wind  tunnel  and  eliminate  the 
negative  aspects  germane  to  the  previous  system.  A  software  package  known  as 
Lab  VIEW©  was  found  to  have  the  attributes  most  needed  to  overcome  the  previous 
system's  inadequacies,  and  had  already  been  purchased  for  another  application.  What 
remained  would  be  to  develop  an  application  that  would  fulfill  the  requirements  of  a  wind 
tunnel  data  acquisition  system. 

The  goals  kept  in  mind  when  developing  a  Lab  VIEW©  application  for  the  wind 
tunnel  were  to  provide  the  user  with  a  well-documented  and  easily  manipulated  program. 
The  program  is  fully  compatible  with  and  functions  under  the  Windows©  operating 
system,  allowing  data  retrieval  onto  portable  media  that  can  be  processed  on  any 
home/office  PC.  Even  more  important,  however,  was  the  need  to  make  results  available 
to  the  user  while  the  data  is  being  acquired.  Real  time  display  of  all  intermediate  steps,  as 
well  as  final  results,  allows  the  user  to  troubleshoot  areas  of  concern  and  pinpoint  sources 
of  technical  problems  regarding  the  testing.  Having  all  intermediate  steps  available  also 
allows  the  user  to  determine  how  various  aspects  of  the  data  contribute  to  the  final  non- 
dimensionalized  results. 

B.  LAB  VIEW©  GRAPHICAL  INSTRUMENTATION  PROGRAM 

Lab  VIEW©  is  a  program  development  application,  much  like  various  commercial 
development  systems  such  as  C  or  BASIC.  However,  LabVIEW©  is  different  from  those 
applications  in  one  important  respect.  Other  programming  systems  use  text-based 
languages  to  create  lines  of  code,  while  LabVIEW©  uses  a  graphical  programming 
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language,  G,  to  create  programs  in  block  diagram  form  [Ref.  1].  The  implication  is  that 
once  designed,  a  Lab  VIEW©  program  will  offer  a  simple  method  for  acquiring  and 
processing  data  obtained  from  the  wind  tunnel  or  other  devices. 

Lab  VIEW©  relies  on  graphical  symbols  rather  than  textual  language  to  describe 
programming  actions.  There  are  extensive  libraries  of  functions  and  subroutines  included 
with  the  software  package  for  many  programming  tasks.  LabVIEW©  contains 
application  specific  libraries  for  data  acquisition,  instrument  control  and  analysis.  Also 
included  in  the  software  are  many  conventional  program  development  tools  to  set 
breakpoints,  animate  program  execution  to  see  how  data  passes  through  the  program,  and 
single-step  through  the  program  to  make  debugging  and  program  development  easier. 

LabVIEW©  programs  are  called  virtual  instruments  (Vis)  because  their 

/  .  . 

appearance  and  operation  imitate  actual  instruments.  Vis  have  both  an  interactive  user 
interface  and  a  source  code  equivalent,  and  accept  parameters  from  higher-level  Vis.  The 
following  are  descriptions  of  these  VI  features  as  taken  from  the  LabVIEW©  User 
Manual  [Ref.  1]. 

1 .  Vis  contain  an  interactive  interface  between  the  user  and  software  which  is 
called  the  front  panel  because  it  simulates  the  panel  of  a  physical  instrument.  The  front 
panel  can  contain  knobs,  push  buttons,  graphs,  and  other  controls  and  indicators.  Data  is 
acquired  by  the  front  panel  via  the  keyboard  and  mouse;  the  results  can  be  viewed  on  the 
computer  screen. 

2.  Vis  receive  instructions  from  a  block  diagram,  which  is  constructed  in 
LabVIEW©’s  programming  language,  "G."  The  block  diagram  supplies  a  pictorial 
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solution  to  a  programming  problem,  and  depicts  graphically  the  written  code  familiar  to 
most  programmers,  e.g.,  “while  loops,”  “for  loops,”  “if/then  cases,”  “formula  nodes,”  etc. 
The  means  whereby  the  front  panel  items  are  wired  to  the  rest  of  the  program  are  also 
displayed.  In  other  words,  the  block  diagram  contains  the  source  code  for  a  given  VI. 

3.  Vis  use  a  hierarchical  and  modular  structure.  They  can  be  developed  and 
used  as  top-level  programs,  or  as  subprograms  within  other  programs  or  subprograms. 
When  a  VI  is  encapsulated  within  another  VI  it  is  called  a  sub  VI.  The  icon  and  connector 
pane  of  a  VI  work  like  a  graphical  parameter  list  so  that  other  Vis  can  pass  data  to  it  as  a 
subVl.  The  above  descriptions  collectively  comprise  what  is  known  as  modular 
programming.  Modularity  in  this  case  leads  to  dividing  an  application  into  a  series  of 
tasks,  which  can  be  further  broken  down  in  a  hierarchical  manner  again  until  a 
complicated  application  comprises  a  number  of  simple  subtasks.  In  the  present  case,  the 
overarching  Vi's  task  is  to  acquire,  process,  display  and  save  force  and  moment  data. 
Various  sub  Vis  acquire  the  voltages,  perform  the  interactions,  plot  the  graphs,  append  to 
spreadsheets,  etc.  Many  of  the  more  complicated  subVIs  contain  subVIs  of  their  own. 
Note  that  it  is  not  mandatory  to  create  subVIs.  The  program  could  execute  as  one  very 
large  and  cumbersome  VI.  However,  since  each  sub  VI  can  be  executed  by  itself  apart 
from  the  rest  of  the  application,  debugging  is  much  easier. 
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C.  TESTING  MOTIVATION 


Two  major  aspects  were  covered  in  this  thesis;  validating  the  LabVIEW© 
program  and  verifying  the  accuracy  of  data  obtained  through  the  LabVIEW© 
methodology  as  compared  to  the  previous  data  acquisition  program.  With  regard  to  the 
LabVIEW©  program,  the  first  step  was  to  confirm  that  the  VI  correctly  yields  the  desired 
results  when  known  values  are  inputted.  A  comparison  was  made  to  the  Quick  Basic 
code  used  by  the  former  data  acquisition  system  to  confirm  similarity  of  results.  Full 
validation  of  the  code  was  accomplished  when  the  VI  demonstrated  the  capability  to 
correctly  acquire  and  process  the  data  from  the  wind  tunnel. 

The  second  aspect  concerns  verifying  the  results  obtained  on  the  Price  Waverider 
by  Cedrun  [Ref.  2].  Based  upon  the  theoretical  potential  of  Price's  design,  NASA  Ames 
constructed  two  aluminum  models  for  testing  Price's  design  at  the  Naval  Postgraduate 
School.  A  15-inch  long  aluminum  model  was  developed  for  wind  tunnel  testing  and  an 
8-inch  long  model  was  made  for  water  tunnel  testing.  Cedrun  tested  the  wind  tunnel 
model  in  1994.  General  results  agreed  favorably  with  HAVOC  predicted  values  for  the 
Price  waverider  and  previous  wind  tunnel  testing  of  Vanhoy's  waverider  and  conventional 
delta  wing  configurations  [Ref.  2].  Questions  did  arise  concerning  pitching  moment 
instability  and  flow  separation,  both  of  which  occurred  at  an  AOA  lower  than  expected. 
In  addition,  axial  forces  were  often  negative,  even  at  low  angles  of  attack.  Therefore  this 
thesis  will  attempt  to  validate  Cedrun's  results  with  the  improved  data  acquisition 
techniques  available  using  LabVIEW©. 
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D.  WHY  LOW-SPEED  TESTING 


The  Naval  Postgraduate  School  (NPS)  Aerolab®  wind  tunnel  is  only  capable  of 
up  to  200  mph  maximum  speed.  In  addition,  the  hardware  used  in  the  tunnel  to  acquire 
data  is  capable  of  handling  only  limited  forces  and  would  be  damaged  by  the  forces 
generated  by  a  higher  speed  tunnel.  Either  of  these  constraints  limits  the  testing 
conditions  to  low  subsonic  flight,  but  this  is  still  sufficient  to  validate  the  performance  of 
the  Lab  VIEW©  VI. 

Waverider  configurations  are  the  result  of  "on-design"  optimization;  that  is,  they 
are  typically  optimized  for  flight  at  specific  design  Mach  number.  Theoretical  flight 
performance  results  are  thus  based  upon  a  specified  set  of  operational  conditions 
including  altitude,  flight  speed  and  a.  Reality  dictates  that  aircraft  flight  is  a  dynamic 
process.  Included  in  this  process  are  maneuvers  or  flight  conditions  that  require  low 
flight  speed,  at  the  very  minimum,  takeoffs  and  landings.  In  general,  there  is  very  little 
.knowledge  of  how  well  suited  the  hypersonic  derived,  geometrical  characteristics  of 
waveriders  are  for  acceptable  subsonic  performance.  This  provides  a  further  rationale  for 
testing  a  hypersonic  configuration  at  lower  subsonic  speeds. 
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II.  EXPERIMENTAL  APPARATUS 


A.  WIND  TUNNEL 

The  NPS  horizontal  low-speed  wind  tunnel  located  in  Halligan  Hall  was  used  for 
conducting  all  experiments.  Manufactured  by  Aerolab®  Development  Company  in  the 
early  1950s,  it  is  a  single  return,  closed  circuit  tunnel  [Ref.  3].  Air  flow  through  the 
tunnel  is  provided  by  a  100-hp  electric  motor  that  drives  a  three-blade,  variable  pitch  fan. 
A  four-gear  transmission  and  a  10:1  contraction  ratio  allow  for  test  section  speeds  of  up 
to  200  miles  per  hour. 

Closed  circuit  wind  tunnels  require  special  considerations  to  reduce  turbulence 
intensity  in  the  test  section.  The  NPS  tunnel  accomplishes  this  through  use  of  stator 
blades  located  directly  behind  the  fan,  two  fine  wire  mesh  screens  six  inches  apart  in  the 
settling  chamber  and  turning  vanes  located  at  each  comer  of  the  tunnel.  As  a  result,  a  test 
section  ambient  turbulence  intensity  of  0.2  %  is  achieved  [Ref.  4], 

The  wind  tunnel  test  section  has  an  area  of  8.75ft2  (45  inches  wide  by  28  inches  in 
height)  and  is  therefore  approximately  l/10th  the  cross  sectional  area  of  the  settling 
chamber.  The  test  section  walls  are  slightly  divergent  fore  to  aft  to  compensate  for  the 
effective  contraction  caused  by  longitudinal  boundary  layer  growth.  Typical  of  most  low 
speed  tunnels,  leakage  occurs  through  the  duct  walls,  which  uncorrected  would  drop 
pressure  below  atmospheric  in  the  test  section.  A  breather  slot  located  immediately 
downstream  of  the  tunnel  test  section  helps  provide  a  return  mechanism  for  lost  air. 
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Adequate  illumination,  visualization  and  access  to  the  test  model  are  provided  by  frosted 
glass  comer  fillet  fluorescent  lights  and  movable  window  tunnel  sidewalls  located  on 
either  side  of  the  test  section.  A  schematic  of  the  tunnel  is  presented  in  Figure  2.1  and  a 
complete  description  of  all  wind  tunnel  characteristics  is  covered  in  the  NPS  Wind 
Tunnel  Lab  Manual  [Ref.  3]. 

Wind  tunnel  air  temperature  is  measured  by  use  of  a  dial  thermometer  extended 

into  the  settling  chamber.  The  test  section  dynamic  pressure,  q  =  i  pxWj,  is.  determined 

from  the  static  pressure  difference,  Ap,  between  four  manifold-flush  static  taps  in  the  test 
section  and  a  similar  set  of  four  taps  in  the  settling  chamber.  Both  sets  of  taps  are 
connected  to  a  common  manifold  and  the  value  for  Ap  is  presented  on  a  micromanometer. 
The  Ap  is  converted  into  dynamic  pressure  using  equation  (2.1)  and  a  calibration  constant 
calculated  from  a  1998  wind  tunnel  calibration  experiment. 

9=  ^P»VW2  =  2.046  KAp  (2.1) 

where:  q  =  dynamic  pressure  (lbf/ft2) 

p„  =  freestream  density  (slugs/fit3) 

Vro  =  freestream  velocity  (ft/s) 

K  =  NPS  wind  tunnel  calibration  constant. 

=  1/0.9012  =  1.109632 

2.046  =  conversion  factor  (cm  H20  to  lbf/ft2) 
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Figure  2.1  NPS  Low  Speed  Wind  Tunnel 
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B.  STRAIN  GAUGE  BALANCE  AND  ASSEMBLY 

Force  and  moment  measurements  were  obtained  by  means  of  a  six  component 
internal  strain  gauge  balance.  The  balance,  on  loan  from  NASA  Ames  and  manufactured 
by  Task  Corporation,  measures  .75  inches  in  diameter  and  is  a  floating-frame  type  as 
pictured  in  Figure  2.2  on  the  following  page.  This  frame  consists  of  an  inner  rod  that  is 
fastened  to  the  support  sting  and  an  outer  casing  that  is  inserted  into  the  waverider  model. 
Forces  and  moments  are  obtained  by  measuring  the  resistance  of  the  elements  connected 
between  the  inner  rod  and  outer  casing. 

All  gauged  sections  of  the  balance  used  during  experimentation  utilize  bending,  4- 
active  foil  type,  epoxy  bonded,  350  ohm  (Q)  strain  gauges.  The  two  normal  force 
elements  (N1  and  N2)  use  350  Q.  bridge  resistance  for  determining  normal  force  and 
pitching  moment,  as  do  the  two -side  force  elements  (SI  and  S2)  for  side  force  and  yawing 
moment.  The  dual  axial  element  (AX)  and  dual  rolling  element  (RM)  each  use  two  full 
wheatstone  bridges  wired  in  parallel  with  a  corresponding  resistance  of  175  Q  for 
determining  drag  and  rolling  moment  respectfully.  The  coordinate  system  within  which 
these  forces  are  defined  is  shown  in  Figure  2.3.  Wiring  for  each  channel  consists  of  four, 
36  gauge  wires,  two  for  the  excitation  voltage  to  the  balance  and  two  for  output  voltage 
from  the  balance  to  the  strain  gauge  signal  conditioner. 
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Figure  2.2  Sting  Balance 


U, 


frees  tream 


Figure  2.3  Sting  Balance  Coordinate  System 
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The  six  output  channels  of  the  strain  gage  balance;  Nl,  N2,  SI,  S2,  AX  and  RM, 
were  calibrated  at  the  NASA  Ames  Calibration  Laboratory  prior  to  release  of  the  balance 
to  the  Naval  Postgraduate  School.  The  excitation  voltage  for  the  calibration  procedure 
was  5  volts  direct  current  (V DC)  with  interaction  constants  determined  iteratively  during 
calibration.  The  interaction  coefficients  determined  by  the  NASA  Ames  calibration  are 
presented  in  Appendix  A.  This  calibration  can  be  verified  by  placing  known  loads  on  the 
balance  and  comparing  them  to  results  displayed  on  the  front  panel  (computer  screen). 

The  balance  is  attached  via  a  6.875  inch  extender  sleeve  to  a  yoke  assembly  that, 
in  turn,  is  mounted  in  the  wind  tunnel  test  section  as  shown  in  Figure  2.4.  Four  machine 
screws  secure  the  bottom  of  the  assembly  to  a  turntable  mechanism  in  the  test  section. 
Depending  on  the  orientation  of  the  model,  the  turntable  facilitates  measurement  of  yaw 
and  angle-of- attack  (AO A).  The  top  of  the  yoke  fits  through  a  circular  cutout  in  the 
ceiling  of  the  test  section,  allowing  rotation  by  the  turntable  but  restricting  any  other 
unwanted  movement.  All  24  wires  from  the  strain  gauge  balance  are  fed  through  the 
sleeve,  yoke  assembly  and  then  out  of  the  tunnel  though  the  breather  (or  ventilation)  slot 
at  the  end  of  the  test  section.  Enough  slack  in  the  wire  bundle  was  allowed  for  ±  90° 
angle  of  rotation. 

* 

C.  STRAIN  GAUGE  SIGNAL  CONDITIONING  BOARD 

Once  outside  the  wind  tunnel,  the  24  wires  from  the  strain  gauge  balance  are 
routed  to  a  National  Instruments  Strain  Gauge  Signal  Conditioning  (SC)  Board,  model  # 
SC-2043-SG.  The  purpose  of  signal  conditioning  is  to  interface  and  properly  prepare  the 


12 


Figure  2.4  Aircraft  Model  Mounting. 


strain  gauge  signals  for  the  data  acquisition  (DAQ)  board  installed  in  the  computer.  The 
SC-2043-SG  provides  for  up  to  eight  channels  of  signal  conditioning  with  a  selection  of 
three  wheatstone  bridge  configurations  (half,  full  and  quarter).  On-Board  excitation  of 
2.5Vdc  is  provided;  however,  it  was  decided  to  use  an  external  excitation  source  of  5VDC 
since  5  was  used  by  NASA  to  determine  the  interaction  coefficients. 

Voltage  excitation  leads  and  input  signal  leads  (i.e.  connections  with  the  strain 
gauge  balance)  are  accomplished  via  a  screw  terminal  block.  The  conditioned  signals 
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from  the  SC-2043-SG  are  transferred  to  the  DAQ  board  via  a  50  pin  connector  and  cable. 
The  cable  connection  also  serves  as  a  power  cable,  providing  power  to  the  SC  board  from 
the  DAQ  board.  A  green  Light  Emitting  Diode  (LED)  indicates  that  the  SC  board  is  on. 

Proper  configuration  of  the  SC  board  is  accomplished  through  one  slide  switch 
and  12  jumpers  (jumper  leads).  The  slide  switch  (SW1)  selects  the  power  supply  for  the 
SC  board,  in  this  setup  the  internal  (INT)  position  is  chosen  to  acquire  power  via  the 
DAQ  board.  Jumpers  JW1  and  JW3  are  used  to  select  onboard  or  external  voltage 
excitation  to  the  strain  gauge  balance.  In  this  case,  the  EXT  position  is  selected. 
Jumpers  JW2  and  JW8  provide  for  the  option  of  excitation  voltage  sensing  on  the  DAQ 
board.  The  remaining  eight  jumpers  configure  the  onboard  (full,  half  or  quarter 
wheatstone)  bridge  completion  network  for  each  of  the  eight  channels  that  may  be  used. 
For  further  details  regarding  configuring  the  SC  Board  consult  the  SC-2043-SG  User 
Manual  [Ref.  5].  Figure  2.5  shows  a  block  diagram  that  illustrates  the  key  functional 
components  of  the  SC-2043-SG. 

D.  AT-MIO-16E-10  DATA  ACQUISITION  (DAQ)  BOARD 

The  AT-MIO-16E-10  is  a  completely  "Plug  and  Play-compatible,"  multifunction 
analog,  digital,  and  timing  I/O  boards  for  use  with  IBM  PCs  and  compatible  computers. 
The  DAQ  board  is  the  interface  between  LabVIEW©  and  the  strain  gauge  signal 
conditioning  board.  It  is  capable  of  acquiring  analog  signals  from  thermocouples,  strain 
gauges,  voltage  sources,  and  current  sources.  It  can  also  acquire  or  generate  digital 
signals  for  communication  and  control.  The  board  includes  a  12-bit  Analog  to  Digital 
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Figure  2.5  SC-2043-SG  Block  Diagram 


5 


50*Pin  I/O  Connector  (or  Lab/1200  Series  Board 


Converter  (ADC)  with  16  analog  inputs,  12-bit  Digital  to  Analog  Converters  (DACs) 
with  voltage  outputs,  eight  lines  of  TTL-compatible  digital  I/O,  and  two  24-bit 
counter/timers  for  timing  I/O  [Ref.  6],  Because  the  DAQ  board  has  no  DIP  switches, 
jumpers,  or  potentiometers,  it  was  easily  configured  and  calibrated  using  the  NI-DAQ 
software. 

A  completely  switchless  and  jumperless  DAQ  board  is  made  possible  by  the 
National  Instruments  DAQ-PnP  bus  interface  chip  which  connects  the  board  to  the  AT 
I/O  bus.  The  DAQ-PnP  implements  the  "Plug  and  Play"  ISA  Specification  so  that  the 
DMA,  interrupts,  and  base  I/O  addresses  are  all  software  configurable.  The  primary 
benefit  of  this  is  that  the  board  configuration  can  be  easily  changed  without  having  to 
remove  the  board  from  the  computer. 

Timing  related  functions  are  handled  by  the  National  Instruments  DAQ-STC 
system  timing  controller.  The  DAQ-STC  consists  of  three  timing  groups  that  include  a 
total  of  seven  24-bit  and  three  16-bit  counters  and  a  maximum  timing  resolution  of  50  ns. 
The  three  groups  control  analog  input,  analog  output,  and  general-purpose  counter/timer 
functions. 
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III.  EXPERIMENTAL  PROCEDURE 

A.  HISTORICAL  PERSPECTIVE 

Since  the  early  1940s,  the  difficulty  in  designing  aircraft  configurations  that 
develop  high  lift-to-drag  (L/D)  ratios  at  high  supersonic  and  hypersonic  speeds  has 
resulted  in  much  research.  Even  in  the  late  1990s,  the  challenge  to  develop  realizable 
hypersonic  aircraft  persists,  with  many  top  designers  and  engineers  pursuing  alternative 
solutions.  One  plausible  concept  that  has  evolved  from  hypersonic  flight  research  is  the 
waverider  configuration.  Though  primarily  based  on  theory,  the  configuration  shows 
high  L/D  potential  and  is  the  focus  of  considerable  aerodynamic  research  in  both  Europe 
and  the  United  States.  In  particular,  ongoing  waverider  research  is  taking  place  at  NASA 
Ames,  NASA  Langley,  University  of  Maryland,  University  of  Oklahoma,  and  the  Naval 
Postgraduate  School.  The  following  paragraphs  give  a  brief  description  of  the  waverider 
concept  and  a  historical  synopsis  detailing  its  evolution.  This  provides  a  background  for 
the  remainder  of  the  paper. 

A  waverider  is  a  supersonic  or  hypersonic  vehicle  that,  at  the  design  point,  has  an 
attached  shock  wave  all  along  its  leading  edge.  Because  of  this,  the  vehicle  appears  to  be 
riding  on  top  of  its  shock  wave  -  hence  the  term  “waverider”  [Ref.  7].  This  is  in  contrast 
to  a  more  conventional  hypersonic  vehicle,  where  the  shock  wave  is  usually  detached 
from  the  leading  edge.  The  aerodynamic  advantage  of  the  waverider,  as  shown  operating 
at  its  design  point  in  Figure  3.1a,  is  “that  the  high  pressure  behind  the  shock  wave  under 
the  vehicle  does  not  “leak”  around  the  leading  edge  to  the  top  surface”  [Ref.  7]. 
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(a)  WAVERIDER 


Figure  3.1  Waverider  vs.  a  generic  hypersonic  configuration. 
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Thus  "the  flow  field  over  the  bottom  surface  is  contained  and  the  high  pressure 
is  preserved”  [Ref.  7],  In  contrast,  for  the  vehicle  shown  in  Figure  3.1b,  there  is 
communication  between  the  flows  over  the  bottom  and  top  surfaces.  In  this  case  the  high 
pressure  air  “leaks”  around  the  leading  edge  and  the  level  of  pressure  on  the  bottom 
surface  is  reduced,  resulting  in  less  lift.  Because  of  this  leakage,  the  conventional  vehicle 
in  Figure  3.1b  must  fly  at  a  larger  AO  A  to  produce  the  same  lift  as  the  waverider  in 
Figure  3.1a.  This  yields  a  secondary  benefit,  which  is  that  because  the  waverider 
generates  the  same  lift  at  a  smaller  a,  drag  is  significantly  reduced  and  the  L/D  for  the 
waverider  is  considerably  higher  than  that  of  a  conventional  supersonic  or  hypersonic 
configuration.  Typically  the  wave  drag  for  the  weak  attached  shock  is  also  lower  than  the 
wave  drag  induced  by  conventional  aircraft  with  a  detached  shock. 

“Another  characteristic  of  waveriders  is  that  they  are  generated  from  known  flow 
fields  established  by  other  basic  shapes  different  from  that  of  the  waverider  itself’  [Ref. 
7],  If  we  restrict  this  consideration  to  a  3D  wing-body  or  lifting  body,  such  as  the  Price 
configuration  that  will  be  introduced  later,  then  the  lower  and  upper  surfaces  of  the 
waverider  are  usually  derived  from  stream  surfaces  in  the  disturbed  flow  generated  by  a 
cone  (or  other  shape).  Even  though  cone-flow  derived  configurations  are  not  conical,  the 
flow  field  in  the  shock  layer  between  the  conical  shock  and  the  non-conical  body  will 
remain  conical.  The  surfaces  of  the  waverider  thus  become  stream  surfaces  that  exists 
behind  a  conical  shock  wave,  stream  surfaces  that  are  generated  by  streamlines  that  begin 
on  the  shock  surface  itself.  Hence,  the  shock  wave  is,  by  definition,  attached  to  the 
leading  edge  of  the  waverider.  Alternatively,  it  can  be  said  that  the  lateral  edges  of  the 
waverider  ride  on  the  captured  shock  wave  [Ref.  8]. 
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Eggers  at  NASA  Ames  initiated  preliminary  research  into  the  waverider  concept 
in  1955.  During  this  time,  Eggers  designed  and  studied  a  flat-top  wing-body  combination 
in  which  the  body  is  situated  entirely  below  the  wing  [Ref.  9].  This  design  was  based  on 
the  principle  that  the  wing  could  be  tapered  to  follow  its  own  shock  wave. 

A  bomber  was  actually  developed  by  North  American  Aviation  that  took 
advantage  of  Egger's  research  (the  B-70  Valkyrie)  although  it  was  not  inversely  derived 
from  a  known  flowfield  [Ref.  10].  In  this  case  the  configuration  allowed  the  shock  wave 
created  by  the  tapering  underbody  of  the  airplane  to  sweep  back  parallel  to  the  leading 
edge  of  the  wing,  just  behind  it  on  the  lower  surface.  A  natural  phenomenon  of  a  shock 
wave  is  that  it  is  a  compression  wave  across  which  a  large  buildup  of  positive  pressure 
occurs.  This  positive  pressure  field  behind  the  shock  wave  was  superimposed  on  the 
underside  of  the  B-70  wing.  The  pressure  field  is  augmented  by  the  positive  pressure  on 
the  tapered  fuselage  and  tends  to  be  contained  by  the  B-70's  folding  wingtips.  The 
pressure  field  thus  created  and  contained  supports  approximately  30%  of  the  weight  of 
the  air  vehicle  at  cruise  conditions.  This  means  the  airplane  can  fly  at  lower  angle  of 
attack  for  a  given  weight,  thereby  decreasing  the  drag  [Ref.  10].  This  indeed  is  very 
similar  to  the  waverider  concept. 

In  1959,  Nonweiler  introduced  what  was  probably  the  first  true  waverider  design. 
The  design  was  based  on  the  idea  of  a  three  dimensional  body  derived  from  the  flow  field 
behind  a  planar  shock  [Ref.  11],  The  resulting  configuration  was  a  delta  wing  planform 
top  view  with  a  caret  (named  after  the  dictionary  symbol  "A")  shaped  cross  section. 
While  flying  at  the  design  Mach  number,  the  shock  would  be  attached  to  the  leading 
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"Shock 

b.  RESULTING  WING  AND  SHOCK 

Figure  3.2  Nonweiler's  Caret  Wing 
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edges  preventing  spanwise  flow  and  spillage  from  the  lower  to  upper  surface  [Ref.  11]. 
Figure  3.2  shows  Nonweiler's  configuration. 

From  the  1960s  through  the  1980s,  waverider  designers  expanded  upon 
Nonweiler's  research  by  exploring  known  flow  fields  generated  by  right  circular  and 
elliptic  cones.  A  conical  flow  waverider  is  shown  in  Figure  3.3  [Ref.  12].  In  the  late 
1980's,  Bowcutt  and  Anderson  continued  development  of  waverider  designs  based  upon 
L/D  optimization  which  included  viscous  effects  [Ref.  13].  Including  viscous  effects  in 
the  design  analysis  was  an  essential  step  for  turning  theory  into  practical  application. 

The  rationale  behind  including  viscous  effects  in  the  optimization  is  based  on  the 
fact  that  previous  configurations  using  only  inviscid  analysis  tended  to  be  less  realistic 
than  those  configurations  based  on  viscous  flows.  The  skin  friction  drag  associated  with 
the  large  surface  area  greatly  reduces  the  predicted  L/D  ratio,  which  led  many  researchers 
to  lose  interest  in  the  waverider  concept  [Ref.  7].  Including  viscous  effects  in  the  L/D 
optimization  allows  for  the  determination  of  a  waverider  with  maximum  L/D  and  takes 
into  account  tradeoffs  between  wave  drag  and  skin  friction  drag  during  the  process. 

Bowcutt  and  Anderson’s  configuration,  obtained  the  lower  surface  of  their  design 
from  stream  surfaces  in  the  disturbed  flow  generated  by  a  cone.  This  is  typical  of 
previous  designs.  The  upper  surface  differed  in  that  it  used  the  stream  surfaces  of  a 
contracting  cylinder  in  its  derivation  [Ref.  13].  In  1988  Vanhoy  contributed  more 
research  to  including  viscous  effects  in  the  optimization  process  of  waverider  design  by 
performing  low  speed  wind  tunnel  comparisons  of  his  own  Mach  6  L/D  optimized 
waverider  with  a  sharp  edged  delta  wing  platform  [Ref.  14].  Boycutt/Anderson  and 
Vanhoy's  optimized  Mach  6  waverider  designs  are  shown  in  Figure  3.4. 
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b.  RESULTING  WING  AND  SHOCK 


Figure  3.3  Conical  Flow  Waverider  Design 
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A.  BOWCUTT  AND  ANDERSON’S  OPTIMIZED  WAVERIDER 


B.  VANHOY'S  MACH  6  OPTIMIZED  WAVERIDER 


Figure  3.4  Mach  6  Optimized  Waveriders 


One  further  step  taken  in  waverider  design  was  the  development  of  a  mission 
specific  waverider  configuration.  In  1993,  LT  David  Price,  USN,  with  assistance  from 
the  NASA  Ames  Research  Center,  completed  the  optimization  and  performance  analysis 
of  a  hypersonic  (Mach  6)  conical-flow  waverider  for  a  deck-launched  intercept  mission 
[Ref.  15].  Using  the  Waverider  Code  and  Hypersonic  Aircraft  Vehicle  Optimization 
Code  (HAVOC)  developed  by  the  Systems  Analysis  Branch  of  NASA  Ames,  and  taking 
into  account  the  practical  considerations  of  optimum  waveriders  discussed  by  Schindel 
[Ref.  16],  a  hydrocarbon-scramjet  powered  waverider  optimized  for  mission  performance 
was  designed  by  Price.  The  configuration  is  unique  in  that  it  not  only  includes  viscous 
effects  but  is  optimized  to  maximize  the  product  of  L/D  and  Isp  [Ref.  15].  Figure  3.5 
shows  Price's  optimum  configuration. 

B.  WAVERIDER  MODEL 

The  Price  Waverider  model  was  manufactured  in  1994  by  the  NASA  Ames 
Machine  Shop  to  accept  the  strain  gauge  (sting)  balance.  7075  aluminum  alloy  was 
chosen  as  the  material  for  the  model  due  to  its  weight,  rigidity  and  ability  to  maintain  an 
edge.  Properties  of  the  aluminum  alloy  are  listed  in  Table  3.1.  To  avoid  any 
misinterpretation  in  reading  the  table,  the  value  given  for  ultimate  strength  equals 
578,000,000  N/m2.  The  sting  balance  is  mounted  via  a  3/4  inch  diameter  hole  drilled  in 
the  base  of  the  model.  It  is  secured  to  the  model  by  a  1/8  inch  set  screw.  This  set  screw 
is  located  along  the  centerline  of  the  model  10.25  inches  aft  of  the  nose  and  corresponds 
to  the  sting  balance  focal  point  when  the  sting  is  mounted  in  the  model.  Model 
parameters  are  listed  in  Table  3.2. 
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Figure  3.5  Price  Mach  6  Waverider 
TABLE  3.1 


7075  ALUMINUM  PROPERTIES 


Property 

density 

Ultimate 

Strength 

Yield 

Strength 

Young's 

Modulus 

Shear 

Modulus 

Units 

BE2S3HB 

Value 

mam 

hhh 

■MMI 

TABLE  3.2 

MODEL  PARAMETERS 


Length 

15  in 

Span 

13.9375  in 

Weight 

7.65445  1b  . 

Planform  Area 

AR 

1.41 

Material 

7075  A1 

C.  DATA  ACQUISITION  PREPARATION 

1.  Strain  Gauge  Signal  Conditioning  Board  Preparation 

With  no  power  to  the  signal  conditioning  board,  a  connection  was  made  between 
it  and  the  DAQ  board  installed  in  the  computer.  A  50  to  22  pin  connector  cable  made  by 
National  Instruments  was  used  to  make  the  connection.  Setup  and  recognition  of  the 
conditioning  board  required  running  the  NI-DAQ  installation  setup  program  included 
with  LabVIEWO.  Under  "configure  accessories"  the  SC-2043-SG  was  selected  with 
voltage  sensing  chosen  as  channel  eight.  (Allowable  strain  gauge  channels  are  numbered 
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0-7,  while  channel  8  is  provided  to  sense  and  display  the  excitation  voltage  on  the  front 
panel  of  the  Lab  VIEW©  VI.)  The  DAQ  board  was  further  configured  as  "nonreferenced 
single-ended"  as  required  by  the  conditioning  board  [Ref.  4]. 

Switch  #1  (SW1)  was  set  to  the  internal  (INT)  position  to  provide  power  to  the 
conditioning  board  from  a  +5VDC  line  on  the  DAQ  board.  This  voltage  source  powers 
the  board  only  and  is  not  used  to  provide  excitation  to  the  strain  gauges.  The  board  does 
provide  a  fixed  excitation  of  2.5 Vdc  but  this  is  insufficient  for  excitation  as  the  sting 
balance  was  calibrated  at  5 Vdc,  therefore  an  external  excitation  source  was  used.  This 
required  setting  jumpers  W1  and  W3  to  the  external  (EXT)  position  and  connecting  the 
excitation  source  to  connector  J7.  Setting  jumper  W8  to  the  EXT  position  allowed  the 
previously  mentioned  voltage  sensing. 

Wheatstone  bridge  completion  was  the  next  item  that  required  configuration, 
although  this  was  simply  verification  since  the  factory  default  settings  were  for  a  full- 
bridge  setup;  the  same  as  the  strain  gauge  balance.  Choices  are  possible  for  full,  half  or 
quarter-bridge  completion,  so  it  was  necessary  to  ensure  that  default  settings  were 
selected.  Jumpers  W4-W7  and  W10-W13  were  verified  to  be  in  the  full-bridge 
configuration. 

The  final  connections  involved  connecting  the  24  wires  from  the  sting  balance  to 
the  screw  terminal  block  mounted  on  the  signal  conditioning  board.  Six  of  the  eight 
channels  of  the  board  were  required  for  the  sting  balance  connections.  Each  channel  of 
the  sting  has  four  wires  each,  two  for  +/-  excitation  (EX)  and  two  for  the  +/-  return  (CH) 
signal.  Since  the  channels  provided  for  this  purpose  are  numbered  0-7  on  the  board  it 
was  arbitrarily  decided  to  use  channels  1  -6  to  keep  things  simple.  They  were  wired  in 
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the  following  order:  Nl,  N2,  SI,  S2,  AX,  and  RM.  Two  channels  on  the  board  remained 
unused.  To  prevent  these  unused  channels  from  saturating,  causing  significant  power 
consumption  and  possible  loss  of  accuracy,  the  screw  terminals  of  the  CH  +/-  inputs  were 
shorted  to  their  respective  negative  excitation  screw  terminal  [Ref.  4], 

The  calibration  of  the  conditioning  board  setup  required  having  the  Lab  VIEW© 
program  activated.  Each  channel  used  by  the  sting  balance  required  nulling,  which 
amounts  to  setting  the  strain  gauge,  return  signal  voltages  to  zero  when  no  load  is  applied 
to  the  balance.  To  provide  for  nulling,  each  channel  on  the  board  has  its  own  circuit  with 
a  trimming  potentiometer  (pot)  to  adjust  the  voltage,  While  viewing  the  voltage  level  on 
the  monitor,  the  pots  were  rotated  until  zero  voltage  on  all  six  channels  was  obtained. 

2.  Lab  VIEW©  VI  Development 

The  most  time  consuming  process  during  the  evolution  of  this  thesis  involved 
creating  the  Virtual  Instrument  for  performing  experimentation.  With  very  little 
Lab  VIEW©  experience  among  the  faculty,  it  was  known  from  the  outset  that  there  would 
be  a  steep  learning  curve  for  the  writer  in  developing  a  program  that  would  not  only 
work,  but  would  be  easily  used  by  the  rest  of  the  department.  To  realize  the  greatest 
benefit  from  using  a  graphical  program,  the  user  should  be  able  to  get  a  "full  solution"  in 
real  time  and  not  have  to  go  "off-line"  to  process  the  data  taken  during  experimentation. 
In  this  particular  example  there  was  a  desire  to  determine,  in  real  time,  the  non- 
dimensional  values  of  Cl,  Cd,  and  Cm  as  a  function  of  the  angle-of-attack  as  the  data  is 
taken.  This  was  achieved,  and  in  addition,  plots  of  Cl  vs.  a,  Cd  vs.  a.  Cm  vs.  a  and  Cl 
vs.  Cd  were  also  rendered  as  the  data  was  taken.  Finally  the  data  is  appended,  via  the 
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LabVIEW©  VI,  to  any  spreadsheet  or  text  program  desired,  e.g..  Excel,  for  future 
analysis.  The  general  procedure  for  developing  this  program  will  be  discussed  in  the 
next  paragraphs. 

The  first  step  in  the  procedure  was  simply  to  acquire  the  voltages  from  the  sting 
balance  and  display  them  on  the  front  panel.  This  was  accomplished  using  a  pre¬ 
designed  VI  from  the  built  in  library.  There  were  several  acceptable  choices  but  the  one 
copied  to  the  block  diagram  was  the  AcquireNscans.vi.  This  VI  automatically 
determines  the  number  of  channels  being  scanned  and  acquires  the  voltages  for  each. 
Front  panel  controls  allow  tailoring  the  length  of  time  the  channels  are  scanned,  the 
number  of  samples  taken  and  the  gain  adjustment  of  the  DAQ  board,  among  others. 
Appendix  A  depicts  the  Front  panel  and  Block  Diagram.  In  this  case,  six  channels  were 
scanned  from  the  sting  balance  and  one  extra  channel  (CHN  8)  for  verifying  the 
excitation  voltage.  Default  values  are  currently  set  at  10,000  scans  over  10  seconds  with 
gain  set  to  1000  via  high  and  low  limit  selections 

Discounting  the  excitation  sensing,  what  were  acquired  were  60,000  voltage 
signals.  Each  of  the  six  channels  is  a  column  array  of  10,000  values.  What  was  needed 
was  a  time  average  of  the  columns,  which  yields  six  voltage  values  for  N1  through  RM. 
This  was  accomplished  by  transposing  the  [10,000  X  6]  matrix  and  then  using  a  subVI  to 
calculate  the  mean  value  of  each  row.  A  "For  Loop"  performs  the  mean  calculation  for 
each  row  sequentially.  The  results  of  this  calculation  are  displayed  on  the  front  panel 
under  the  title  "Mean."  When  no  load  is  applied  to  the  sting  balance  the  mean  values 
should  read  zero;  therefore,  nulling  the  system  requires  monitoring  these  values  as 
discussed  in  the  previous  section. 
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The  next  step  in  the  block  diagram  requires  a  calculation  of  the  direct  force 
nonlinear  values  and  the  force  interaction  values.  Two  separate  subVIs  were  created  to 
accomplish  this,  preinter  act. vi  and  interact  l.vi.  For  the  first  sub  VI,  the  six  mean  values 
described  in  the  paragraph  above  were  sent  to  individual  formula  nodes  that  execute 
calculations  for  both  positive  and  negative  voltages.  A  "True/False"  case  structure  is 
"wired"  in  the  sub VI  to  the  incoming  signal  from  each  channel  to  determine  if  it  is 
positive  or  negative,  which  determines  which  equation  to  calculate.  This  yields  12 
equations,  two  per  channel.  Appendix  B  lists. the  constants  for  the  direct  force  nonlinear 
equations. 

The  calculations  from  the  preinter  act. vi  are  then  sent  to  inter  act  l.vi  to  account  for 
force  interactions  between  the  strain  gauge  elements  in  the  sting  balance.  This  VI 
executes  calculations  in  much  the  same  way,  with  a  formula  node  for  each  channel  again 
surrounded  by  a  "True/False"  case  structure  to  account  for  positive  or  negative  voltage 
cases.  This  again  yields  12  equations.  However,  since  the  VI  is  checking  for 
convergence  due  to  interactions,  each  calculated  signal  is  routed  back  as  an  input  to  the 
five  other  channels.  Including  the  initial  signal  from  the  preinter  act. vi  there  are  six 
voltage  signals  routed  in  and  one  voltage  signal  routed  out  for  each  channel.  A  For  Loop 
surrounds  the  interact  l.vi  to  perform  the  iterations  with  the  default  number  of  iterations 
set  to  20  on  the  front  panel.  This  seems  enough  to  ensure  convergence  to  the  actual 
value,  as  iterations  greater  than  10  effected  no  change  on  the  calculated  values  during 
testing. 

The  six  resulting  values  from  the  interaction  convergence  yield  the  actual  forces 
felt  by  the  strain  gauge  elements.  These  forces  are  displayed  on  the  front  panel  as  XN1, 
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XN2,  XS1,  XS2,  XAX  and  XRM.  The  next  step  is  to  convert  these  individual  element 
forces  into  forces  that  the  model  "feels."  This  is  accomplished  by  yet  another  formula 
node  in  the  block  diagram.  It  is  also  at  this  point  that  one  final  "True/False"  case  is 
needed.  The  "true"  case  will  select  the  "Tare  reading,"  which  is  simply  the  weight  of  the 
model  with  no  external  forces  applied.  The  "false"  case  will  select  the  forces  on  the 
model  when  the  wind  tunnel  is  running  and  execute  the  remainder  of  the  program. 

The  Tare  reading  is  selected  prior  to  running  the  tunnel.  This  is  achieved  simply 
enough  by  selecting  "ON"  for  the  "TARE"  button  on  the  front  panel.  The  formula  node 
then  computes  the  normal,  side,  axial  and  rolling  (moment)  forces  on  the  model.  These 
forces  should  be  negligible  except  for  the  side  force.  (Recall  that  the  model  is  mounted 
sideways  as  shown  in  Figure  2.4  and  Appendix  E.)  The  calculated  values  are  saved  as 
.local  variables  so  they  will  be  remembered  during  wind  tunnel  testing.  Additionally, 
these  values  are  displayed  on  the  front  panel  as  "Tare  Readings". 

After  the  tunnel  is  running  and  stabilized,  data  collection  is  initiated  by  selecting 
OFF  for  the  Tare  reading.  This  changes  the  "True/False"  case  allowing  force  calculations 
to  be  made.  Other  required  selections  at  this  time  include  entering  AO  A,  manometer  Ap, 
chord  length,  temperature,  atmospheric  pressure  and  planform  area,  as  well  as  selecting  a 
file  to  save  the  data  to  for  future  retrieval.  Selecting  the  run  program  allows  the  formula 
node  to  execute  again,  this  time  calling  the  local  variables  from  the  Tare  readings  and 
subtracting  them  from  the  force  calculations.  The  adjusted  force  calculations  are  then 
collected  as  an  array  and  displayed  on  the  front  panel  below  the  Tare  readings.  They 
consist  of  normal,  axial  and  side  force,  as  well  as  pitch,  yaw  and  rolling  moments. 
Monitoring  the  force  calculations  on  the  front  panel  serves  as  a  useful  tool  in  avoiding 
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overloading  the  balance.  Warning  indicators  turn  red  on  the  front  panel  when  forces 
exceed  80%  of  rated  strain  gauge  loads  to  further  reduce  the  possibility  of  overload. 

As  stand  alone  values,  the  force  calculation  values  are  not  too  useful  to  the 
aerodynamicist  until  they  are  nondimensionalized.  The  next  two  formula  nodes 
accomplish  this;  the  first  calculating  CN,  CA  and  CM  from  the  normal  force,  axial  force 
and  pitching  moment.  Included  in  the  first  formula  node  are  the  user's  inputs  for  chord 
length  and  planform  area,  as  well  as  the  value  for  dynamic  pressure  calculated  from  the 
manometer  Ap.  The  second  formula  node  uses  the  values  of  Cn  and  CA,  along  with  the 
user's  input  for  AO  A,  to  calculate  Cl  and  Cd- 

The  calculated  values  from  the  two  formula  nodes  used  by  this  program  are  AO  A, 
Cl,  Cd  and  Cm-  These  outputs  are  then  collected  into  an  array  in  the  order  stated  above 
and  sent  to  another  sub VI  that  appends  the  data  by  row  to  a  spreadsheet  file.  Other  data 
appended  automatically  to  the  spreadsheet  includes  Cn,  Ca,  Normal,  Axial  and  Side 
forces,  Pitch,  Yaw  and  Rolling  moments,  IAS,  True  wind  tunnel  speed,  Reynolds  number 
and  dynamic  pressure.  Each  successive  run  adds  another  row  to  the  spreadsheet.  It  is 
worth  noting  here  that  other  values  based  on  the  force/moment  calculations  can'  be 
calculated  and  displayed  with  little  effort  if  desired. 

Finally,  the  results  are  graphically  displayed  on  the  front  panel  as  plots  of  Cl  vs. 
a,  Cd  vs.  a,  Cm  vs.  a  and  Cl  vs.  Cd-  A  buffer  sub VI  is  placed  in  front  of  each  graph 
icon  on  the  block  diagram  to  retain  each  data  point  taken.  Otherwise  each  successive  run 
would  overwrite  the  previous  data  point.  The  benefit  of  creating  plots  in  real  time,  in 
addition  to  saving  the  data,  is  a  major  advantage  of  the  program.  The  user  can  monitor 

the  plots  and  verify  that  the  model  is  behaving  as  predicted.  Unexpected  results  can  be 
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immediately  determined,  and  if  possible,  corrected  on  the  spot,  eliminating  repeat  tests 
due  to  faulty  data  collection  or  procedures. 

D.  STING  BALANCE  WIND  TUNNEL  EXPERIMENTS 

1.  Test  plan 

Following  calibration  of  the  signal  conditioning  board  and  sting  balance  [App.  A], 
the  model  was  secured  on  the  yoke  in  the  wind  tunnel.  The  determined  schedule  of 
testing  was  six  data  runs  at  different  tunnel  speeds,  varying  the  AOA  during  each  run. 
The  sweep  ranges  of  a  were  dictated  by  the  rated  load  of  the  balance.  To  ensure  a 
margin  of  safety,  applied  loads  were  not  allowed  to  exceed  80%  of  the  rating.  To  further 
reduce  the  possibility  of  damage  to  the  balance,  each  test  run  was  begun  at  the  lowest 
practical  tunnel  speed.  Tunnel  speed  was  increased  gradually  for  each  run,  with  the  force 
equations  continuously  monitored  for  near  overload  conditions.  A  synopsis  of  the  runs 
is  listed  in  Table  3.1. 

2.  Tare  Values 

Prior  to  taking  wind  tunnel  measurements,  the  tare  values  representing  forces  and 
moments  due  to  the  weight  of  the  model  were  determined.  A  discussion  of  how  this  is 
accomplished  is  covered  under  LabVIEW©  VI  development.  In  short  this  simply  means 
selecting  the  "ON"  button  for  the  TARE  reading.  Subsequent  runs  will  automatically 
subtract  this  amount  from  the  "tunnel  on"  data  points.  It  was  found  that  the  nulling  range 
of  the  circuit  board  allowed  nulling  the  Tare  values  down  to  zero. 
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TABLE  3.3 


TEST  CONDITIONS 


Test 

q  (lb/ft2) 

IAS  (mph) 

Sweep  range  (deg) 

0 

3.61 

45 

0  to  45 

1 

3.65 

45 

-90  to  90 

2 

7.74 

-25  to  22 

3 

13.74 

75 

-18  to  13 

4 

20.29 

-13  to  8 

5 

25.92 

100 

-11  to  6 

6 

31.94 

110 

-10  to  4 

3.  Preliminary  Test 

A  preliminary  test  (test  0)  was  performed  neglecting  the  use  of  the  interaction 
equations.  This  was  done  to  isolate  any  problems  with  data  acquisition  using  the 
Lab  VIEW©  program.  Preliminary  procedures  before  beginning  the  test  were  performed 
to  ensure  integrity  of  the  tunnel  [Ref.  3].  Tunnel  speed  was  brought  up  to  45  mph,  with  a 
corresponding  dynamic  pressure  of  3.61  lbf/ft2.  Angle-of-attack  was  varied  from  0-45° 
while  maintaining  a  constant  vigil  on  the  monitor  to  check  balance  loads.  All  loads  were 
significantly  below  warning  levels,  indicating  much  higher  speed  was  possible.  Data 
acquisition  performed  as  expected,  with  both  real  time  plots  and  spreadsheet  data  created. 
Negative  sweeps  were  to  be  performed  but  the  turntable  drive  mechanism  that  controls 
changing  AOA  failed.  This  temporarily  halted  further  testing  until  repairs  could  be  made 
on  29  August  1998. 
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4.  Subsequent  Tests: 

Tests  one  through  six  comprised  the  formal  testing  of  the  model  with  force 
interaction  coefficients  included.  For  each  run,  the  model  was  initially  placed  at  0°  a 
then  moved  incrementally  to  locate  the  angle-of-attack  that  yielded  forces  equaling  or 
exceeding  the  80%  limit  load  criteria.  It  was  found  in  all  cases  that  N1  reached  its  load 
limit  first,  and  higher  negative  a  settings  were  achievable  than  positive.  Once  the 
limiting  negative  AOA  was  determined  the  model  was  moved  incrementally  (via  the 
turntable)  through  the  allowable  angle  of  attack  range.  Data  points  were  taken  statically 
after  moving  the  turntable  to  each  desired  AOA  location. 

5.  Data  reduction 

Data  files  of  all  tests  were  collected,  as  previously  mentioned  in  spreadsheet 
format,  and  are  presented  in  Appendix  D.  The  LabVIEW©  program  calculated  values  of 
the  force  coefficients  from  the  dimensional  forces  in  the  following  manner.  The  normal 
and  axial  force  coefficients  (Cn  and  Ca)  were  calculated  using  the  appropriate  form  of 
equation  (3.1). 


Where:  Cf  -  force  coefficient 

F  =  force  (lbf) 
q  =  dynamic  pressure  (lbf/ft2) 
S  =  planform  area  (ft2) 


(3.1) 
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Using  equations  (3.2)  and  (3.3),  the  values  for  CN  and  CA  were  inputted  to  calculate  the 
coefficients  for  lift  and  drag.  Figure  3.6  depicts  the  use  of  these  equations. 


Cl=Cn cos  a- CA sin  a 

(3.2) 

Cd  =  Cn  sin  a  +  CA  cos  a 

(3.3) 

Where: 

Cn  =  normal  force  coefficient 

CA  =  axial  force  coefficient 

a  =  angle  of  attack  (degrees) 

Figure  3.6  Calculation  of  Cl  and  Cd 
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The  pitching  moment  coefficient  was  calculated  using  equation  (3.4). 


Cm — 


M 
qS  c 


Where:  Cm  =  moment  coefficient 

M  =  moment  (ff-lbf) 
c  =  root  aerodynamic  chord  (ft) 


(3.4) 


The  length  of  the  model,  i.e.  root  chord  of  1.25  ft,  was  used  as  the  reference  chord  length. 
Plots  obtained  from  these  calculations  include  CL,  CDand  CM  vs.  a,  CLvs.  CD>  CMvsCL 
and  Gl/Cd  vs  a.  The  results  of  these  plots  are  summarized  in  the  next  chapter. 
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IV.  RESULTS  AND  DISCUSSION 


A.  VALIDATION  OF  LAB  VIEW©  PROGRAM 

Before  obtaining  results  from  the  wind  tunnel  testing,  it  was  necessary  to  validate 
the  code  that  was  written  to  acquire  and  process  the  data.  This  was  accomplished  by  first 
running  the  program  offline  and  inserting  known  values  in  the  calculations,  then  by 
acquiring  data  via  the  strain  gauge  balance  and  comparing  results  to  expected  values. 

1.  Offline  Validation 

A  version  of  the  Quick  Basic  code  used  by  Cedrun  [Ref.  2]  was  obtained  and 
converted  to  PC  format,  allowing  manual  input  of  dummy  voltage  values  to  be  entered. 
Dummy  voltages  for  the  six  strain  gauge  channels  (Nl-RM)  were  entered  yielding 
Normal,  Side  and  Axial  “force”  results  and  Pitch,  Yaw  and  Rolling  “moment”  results. 
These  same  dummy  voltages  were  inputted  into  the  Lab  VIEW©  program  and  identical 
results  were  obtained  when  using  Cedrun’ s  interaction  coefficients. 

2.  Acquiring  Data 

With  setup  accomplished  in  the  wind  tunnel,  known  loads  were  applied  in  weight 
increments  of  one  to  ten  pounds  on  the  strain  gauge  balance  as  discussed  in  Appendix  B. 
Results  displayed  on  the  front  panel  confirmed  that  accurate  data  was  obtained  for  the 
weights  applied.  This  served  a  secondary  purpose  of  verifying  the  NASA  Ames 
calibration  of  the  balance. 
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B.  WIND  TUNNEL  DATA 


Results  of  the  wind  tunnel  testing  are  decomposed  into  two  areas.  Lift  and  Drag 
comprising  the  first  part  of  the  discussion  and  pitching  moment  the  second  part.  In  both 
cases  tunnel  wall  corrections  were  applied  to  the  data.  These  corrections  were  made  to 
account  for  the  change  in  flow  properties  caused  by  the  walls  of  the  wind  tunnel's  closed 
test  section.  Corrections  were  included  for  solid  blocking,  wake  blocking  and  streamline 
curvature.  A  brief  discussion  regarding  these  corrections  and  why  they  were  needed  is 
included  in  Appendix  D.  For  more  information,  [Ref.  21]  provides  a  detailed  analysis  of 
wind  tunnel  wall  corrections  and  the  theory  behind  them. 

1.  Lift  and  Drag 

The  graphs  of  Cl  versus  angle  of  attack  for  all  six  test  runs  are  shown  in  Figures 
4.1  through  4.6.  The  graphs  show  little  variation  in  Cl  with  respect  to  dynamic  pressure. 
The  graph  for  Cl  in  Figure  1  is  worth  special  note  due  to  the  large  amount  of  difference 
between  uncorrected  values  and  corrected  after  20°  AOA.  Applying  the  tunnel  wall 
corrections  yields  a  significant  loss  of  CL  for  AOAs  greater  than  20°.  Tunnel  wall 
corrections  for  values  under  20°  AOA  were  minor  in  comparison,  which  was  also  true  for 
all  other  coefficients  plotted.  A  maximum  loss  of  0.5  (CL  drops  from  1.5  to  1.0)  occurs  at 
approximately  30°  AOA,  at  which  point  the  waverider  stalls.  Figure  4.1  substantiates 
Cedrun's  conclusion  that  the  waverider  continues  to  generate  lift  at  higher  angles  of 
attack  after  stall  [Ref.  14]. 

Figures  4.7  through  4.12  display  graphs  for  CD  versus  AOA.  It  can  be  seen  that 
minimum  drag  (CD  =  0.03)  occurs  at  0°  AOA,  therefore  CDmin  =  CDo-  Drag  is  somewhat 
higher  than  results  previously  obtained  by  Cedrun,  which  is  later  shown  to  be  a  result  of 
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the  corrections  made  to  the  interaction  coefficients.  These  corrections  led  to  a  larger 
positive  axial  force  for  all  results  that  in  turn  led  to  higher  drag  values.  The  higher  drag 
characteristics  at  positive  vice  negative  AOA  do  agree  with  Cedrun  and  are  explained  by 
the  corresponding  higher  lift  characteristics  also  observed  at  positive  AOA.  CJ  Cd 
versus  AOA  plots  are  shown  in  Figures  4.13  through  4.18.  As  can  be  seen,  L/D  is  not 
much  affected  by  tunnel  wall  corrections.  A  maximum  value  of  approximately  (L/D)  =  4 
occurred  at  6°  AOA.  Cl  versus  Cd  plots  are  shown  in  Figures  4.25  through  4.30.  Again, 
it  can  be  seen  that  applying  tunnel  wall  corrections  has  a  significant  effect  at  higher 
angles  of  attack,  as  shown  by  Figure  4.25.  From  0  to  +/-  20°,  the  differences  in  Cl  versus 
Cd  are  much  less  apparent. 

A  graph  of  Cl  versus  AOA  attack  with  Reynolds  number  as  a  parameter  is 
provided  in  Figure  4.37.  The  graph  clearly  shows  that  Cl  is  virtually  unaffected  by  a 
change  in  the  Reynolds  number,  also  agreeing  with  previous  results  obtained  by  Cedrun. 

2.  Pitch 

The  graphs  of  Cm  versus  AOA  for  all  six  tests  are  shown  in  Figures  4.19  through 

4.24.  As  with  Cl,  there  was  minimal  variation  in  pitching  moment  with  respect  to 

changing  dynamic  pressure  or  Reynolds  number.  The  point  at  which  the  pitching 

moments  were  initially  taken  is  located  10.25  inches  from  the  nose  on  centerline.  Prior  to 

* 

applying  wind  tunnel  boundary  corrections,  a  moment  transfer  to  the  1/4  root  chord  was 
necessary  due  to  formula  requirements  [Ref.  21].  This  involved  a  transfer  of  6.5  inches, 
locating  the  1/4  chord  at  3.75  inches  aft  of  the  nose  (0.3125  ft). 

Additional  plots  of  Cm  versus  AOA  were  created  to  show  pitching  moment  about 
the  leading  edge  in  Figures  4.31  through  4.36.  This  location  was  chosen  primarily  to 
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compare  with  results  contained  in  an  AIAA  paper  by  Miller  and  Argrow  on  a  high-speed 
civil  transport  (HSCT)  waverider,  which  was  compared  to  the  Price  Waverider  among 
others  [Ref.22],  It  is  suggested  that  the  authors  did  not  transfer  Cedrun's  moment 
location  [Ref.  2],  which  was  at  68%  aft  of  the  LE  of  the  root  chord,  before  making  the 
comparison.  Figure  4.41  displays  a  correct  comparison  of  the  Price  and  HSCT  waverider 
for  Cm  versus  a. 

In  any  case,  the  Price  waverider  appears  longitudinally  stable  when  referenced  to 
the  quarter  chord  and  leading  edge,  but  neither  are  particularly  useful  values  for 
determining  the  longitudinal  stability  of  an  aircraft.  In  order  to  determine  longitudinal 
stability  the  aerodynamic  center  (or  neutral  point),  and  the  aircraft's  center  of  gravity 
(c.g.)  must  be  located.  Since  this  is  a  model,  the  location  of  the  c.g.  is  not  likely  to 
represent  a  real  aircraft,  but  the  aerodynamic  center  of  the  Price  waverider  model  was 
located  at  4.8  in  (32%)  aft  of  the  leading  edge.  For  static  stability,  this  will  require  a  c.g. 
location  less  than  32%  aft  the  leading  edge  of  the  model;  hence  it  appears  unlikely  that  a 
representative  aircraft  will  be  statically  stable.  Figure  4.38  displays  a  plot  of  CM  versus 
AOA  with  zero  slope,  thereby  indicating  that  the  neutral  point  is  at  approximately  CM(C/3). 

The  location  of  the  neutral  point  can  also  be  verified  using  the  following  formula 
if  the  location  of  the  mean  aerodynamic  chord  is  known. 


(4-1) 


The  mean  aerodynamic  chord  was  found  by  determining  the  equation  of  the  line  fitting 
the  perimeter  of  the  waverider  model  as  shown  in  Figure  4.39,  then  applying  equation 
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4.2.  This  defined  the  length  of  the  mean  aerodynamic  chord  to  be  12.05  inches. 


With  reference  to  the  MAC  of  12.05  inches,  its  origin  is  2.95  inches  aft  of  the 
centerline  leading  edge  and  the  1/4  chord  is  located  5.96  inches  from  the  centerline 
leading  edge.  Taking  Cm  versus  Cl  about  the  MAC  1/4  chord  location  yields  the  result 
seen  in  Figure  4.40  when  applying  the  data  from  the  IAS  of  75  mph.  The  slope  of  a 

linear  fit  to  the  Cm  vs  CL  plot  is  0.0978,  which  yields  *  =  0.25-0.0978=  0.1522.  To 

c 

find  the  location  of  the  neutral  point  with  respect  to  the  MAC,  0.1522  was  multiplied  by 

c  yielding  ,  an  x  location  1.83  inches  aft  of  the  leading  edge  of  the  mean  aerodynamic 
chord.  When  referenced  back  to  the  root  chord,  1.83  +  2.95  =  4.78  inches  from  the 
leading  edge  on  centerline  as  can  be  seen  in  Figure  4.39.  This  closely  agrees  with  the  4.8 
inches  found  from  the  previous  method. 

C.  ACCURACY  OF  DATA 

Figures  4.42  through  4.46  compare  data  obtained  from  Cedrun's  work  [Ref.  2] 
with  current  results.  It  was  previously  noted  that  Cedrun  obtained  lower  drag  and 
therefore  higher  L/D  values,  believed  due  to  corrections  made  to  the  interaction 
coefficients  during  the  current  testing.  In  addition,  for  reasons  perhaps  also  due  to  the 
incorrect  interaction  coefficients,  many  of  the  axial  force  values  obtained  by  Cedrun  were 
slightly  negative,  which  would  also  lead  to  lower  drag  values.  However,  his  method  of 
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converting  dynamic  pressure  from  cm  H2O  to  lbf/ft2  yielded  values  of  q  that  were  3.1% 
lower  than  current  values.  This  q  discrepancy  would  tend  to  increase  values  for  both  CD 
and  CL.  With  these  many  differences,  it  was  decided  to  perform  an  additional  test  in  the 
wind  tunnel  in  an  attempt  to  exactly  duplicate  Cednm's  results. 

Interaction  coefficients  identical  to  Cedrun's  were  used  and  testing  was  performed 
using  a  q  obtained  by  his  method.  Results  of  this  final  test  match  almost  exactly  with 
results  Cedrun  obtained.  Closeness  of  fit  is  particularly  apparent  in  Figures  4.43  and  4.45. 
The  curve  fit  of  Cedrun's  data  in  Figure  4.46  is  somewhat  deceiving  as  it  does  not  fall 
exactly  on  line  with  the  LabVIEW©  result.  This  is  due  to  the  limited  number  of  data 
points  for  Cedrun's  test  run.  In  any  case,  this  latest  test  supports  the  validity  of  current 
values  obtained  during  runs  one  through  six.  The  differences  between  these  current 
values  and  those  obtained  during  the  1994  test  data  are  clearly  due  to  corrections  made  to 
the  interaction  coefficients  and  the  method  previously  used  by  Cedrun  in  obtaining  q. 
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Figure  4.1  CL  vs.  AO  A,  IAS  =45  mph 


Figure  4.2  CL  vs.  AOA,  IAS  =  60  mph 
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Figure  4.6  Cl  vs.  AO  A,  IAS  =  110  mph 
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Figure  4.7  CD  vs.  AO  A,  IAS  =  45  mph 


Figure  4.8  CD  vs.  AO  A,  IAS  =  60  mph 
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Figure  4.9  Cd  vs.  AO  A,  IAS  =  75  mph 


Figure  4.10  Cd  vs.  AO  A,  IAS  =  90  mph 
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Figure  4.11  CD  vs.  AOA,  IAS  =  100  mph 


Figure  4.12  Cd  vs.  AOA,  IAS  =  110  mph 
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CL/CD  vs  AOA 


AOA 


Figure  4.13  Cl/Cd  vs.  AOA,  IAS  =  45  mph 

CL/CD  vs  AOA 


AOA 


Figure  4.14  Cl/Cd  vs.  AOA,  IAS  =  60  mph 
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c,/c 


.  Figure  4.15  Cl/Cd  vs.  AO  A,  IAS  —  75  mph 
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C,/Cn  ^  I  C./C 


Figure  4.17  CL/CD  vs.  AOA,  IAS  =  100  mph 


Cl/Cdvs.AOA 


_*_CUCD  vsAOA 
!  _ CL/CD  vs  AOA  corrected 


Figure  4.18  Cl/Cd  vs.  AOA,  IAS  =  110  mph 


Figure  4.19  CM  vs.  AOA,  IAS  =  45  mph 
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Cm(1/4c)  vs  AOA 


Figure  4.20  CM  vs.  AOA,  IAS  =  60  mph 
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Figure  4.21  Cm  vs.  AO  A,  IAS  =  75  mph 
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Figure  4.22  Cm  vs.  AOA,  IAS  =  90  mph 
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Figure  4.23  Cm  vs.  AO  A,  IAS  =  100  mph 


Figure  4.24  CM  vs.  AO  A,  IAS  =  110  mph 


Figure  4.26  Cl  vs.  Cd,  IAS  =  60  mph 
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Figure  4.27  CL  vs.  CD,  IAS  =  75  mph 


Figure  4.28  CL  vs.  CD,  IAS  =  90  mph 


Figure  4.32  Cm(le>  vs.  AOA,  IAS  =  60  mph 
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Figure  4.33  Cm(le)  vs.  AO  A,  IAS  =  75  mpb 


Figure  4.34  Cm(le>  vs.  AO  A,  IAS  =  90  mph 
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Figure  4.35  CM(le)  vs.  AO  A,  IAS  =  100  mph 


Figure  4.36  CM(le>  vs.  AO  A,  IAS  =  110  mph 


CL  vs.  AOA  (Re  independence) 


Figure  4.37  CL  vs.  AOA  (Reynolds  Independence) 


Figure  4.38  Neutral  point  at  1/3  of  root  chord 
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Length 


Waverider  Geometry  and  MAC  location 


MAC  =  2*  integral  of  c2  (semispan) 
planform  area 

=  1721.62  in  /143.35  in  =12.01  in 

MAC  =  1  ft 


- - -6-^- 


Planform  Area  =  .99548  ft2 


jc  =  -0.0007X4  +  7E-15x3  -  0.2653x2  -  2E-12x  +  14.909j 
1  R2  =  0.9992  j 


Figure  4.39  Planform  Geometry,  Location  of  Mean  Aerodynamic  Chord 


Figure  4.40  Determining  slope  of  CM(maq  vs.  Cl 


Figure  4.42  Cl  vs.  AOA  comparison,  Huff  vs.  Cedrun 
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Figure  4.43  Cd  vs.  AOA  comparison,  Huff  vs.  Cedrun 


Figure  4.44  Cm  vs.  AOA  comparison,  Huff  vs.  Cedrun 
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Figure  4.45  Cl  vs.  Cd,  Huff  vs.  Cedrun 


V.  CONCLUSION  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  use  of  a  modem  data  acquisition  system  such  as  Lab  VIEW©  has  significantly 
improved  the  capability  of  the  NPS  Low  Speed  Wind  Tunnel.  When  combined  with  the 
strain  gauge  balance,  the  aerodynamic  forces  on  any  sting-mounted  model  can  be  quickly 
and  accurately  determined.  Several  additional  benefits  are  realized.  First,  data  is 
acquired  and  processed  online  so  results  can  be  analyzed  while  testing  is  taking  place. 
Second,  modifications  to  the  program  can  be  made  easily,  allowing  a  greater  variety  of 
testing  to  be  accomplished  in  a  shorter  amount  of  time.  Third,  the  user  can  monitor  the 
data  processing  on  the  front  panel  through  all  the  intermediate  steps  as  well  as  simply 
viewing  the  final  result.  Finally,  the  program’s  scan  rate  can  be  adjusted  from  the  front 
panel  instantly,  allowing  the  user  to  choose  between  static  testing  for  single  run  program 
execution  or  dynamic  testing  on  a  continuous  running  cycle. 

The  use  of  a  PC  based  system  has  other  inherent  advantages.  Post  processing,  if 
desired  for  future  use,  is  made  easier  by  compatibility  with  existing  spreadsheet  and  word 
processing  programs.  The  data  can  be  easily  archived  on  the  hard  drive  and  accessed 
through  Windows©  Explorer.  Processed  data,  including  all  intermediate  calculations  as 
well  as  final  results,  can  be  saved  automatically  to  floppy,  zip  or  any  other  portable  media 
with  later  retrieval  on  any  other  PC.  The  entire  strain  gauge  program,  due  to  its 
modularity  and  method  of  calling  Vis  from  the  built  in  software  library,  is  also  easily 
portable  and  can  be  saved  on  a  single  floppy  diskette. 

Testing  performed  on  the  Price  Waverider  model  accomplished  several  purposes. 

Data  obtained  during  current  testing  validated  previous  research  performed  by  LT 
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Cedrun.  Current  testing  also  resolved  questions  regarding  longitudinal  stability  of  the 

Price  waverider  and  the  location  of  its  neutral  point.  Perhaps  most  important,  using  a 

model  that  had  been  previously  tested  served  as  a  verification  process  for  the  Lab  VIEW© 

program,  paving  the  way  for  future  wind  tunnel  analysis  of  other  models. 

B.  RECOMMENDATIONS 

1 .  The  Price  waverider  has  been  evaluated  so  far  with  an  emphasis  on  static  testing.  The 
capability  exists  for  dynamic  and  oscillatory  testing  as  well. 

2.  The  range  of  AO  A  was  quite  limited  above  100  mph  IAS  using  the  larger  wind  tunnel 
model.  A  water  tunnel  model  of  the  Price  waverider  can  be  tested  in  the  wind  tunnel 
using  the  locally  manufactured  sting  adapter.  As  the  planform  area  is  significantly 
smaller,  testing  can  be  conducted  at  higher  speeds  without  fear  of  overloading  the 
balance.  A  comparison  could  be  made  (in  the  wind  tunnel)  between  the  wind  tunnel 
and  water  tunnel  models  to  verify  similarity  of  results. 

3.  Other  water  tunnel  models,  including  another  waverider,  are  also  available  for  testing 
in  the  wind  tunnel  using  the  same  locally  manufactured  adapter.  Comparisons  of 
aerodynamic  qualities  could  be  easily  made  between  these  different  models. 

4.  The  wind  tunnel  is  capable  of  200  mph,  far  exceeding  the  capabilities  of  the  current 
strain  gauge  balance  with  larger  models.  The  use  of  a  balance  with  higher  load 
ratings  would  allow  larger  changes  in  AOA  at  increased  tunnel  speeds. 

5.  Side  force,  Yaw  and  Rolling  moments  were  not  non-dimensionalized.  Future 
modifications  could  easily  incorporate  these  features  into  the  existing  program. 
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APPENDIX  A:  LABVIEW©  PROGRAM 


The  first  three  pages  of  Appendix  C  comprise  the  front  panel  of  the  program. 
This  is  where  the  user  makes  inputs  and  views  the  results.  The  remainder  of  the  appendix 
displays  the  overall  block  diagram  and  block  diagrams  of  the  preinteract  and  interactl 
sub  Vis.  The  block  diagrams  show  the  process  from  acquisition  to  final  analysis/plotting. 


General  Guidelines: 

L  TAPE  Calculation 

1.  Select  "OFF"  for  append  to  existing  file. 

2.  Type  in  a  file  destination  for  the  data. 
This  data  is  fore  spreadsheet,  ie.  Excel. 

Each  row  will  yield  one  value  of  AO  A, 

Cl  Cd  end  Cm  (Suggest  usea  floppy,  in  which 
case  you  could  type  A  \data.xls) 

3.  Leave  "format”  and  "delimiter"  alone. 


4.  Select  "ON"  for  T are  reading. 

5.  Choose  #  of  Iterations  for  the 
Interaction  loop".  10  is  a  g ood  number. 

S.  Select  "ON"  for  real  data  collection. 

7.  Take  Tare  calculation  by  clicking  arrow 
in  upper  left  comer  of  toolbar. 

1  Remaining  Calculations 

1.  Select  "ON"  for  append  to  existing  file. 

2.  Select  “OFF"  for  Tare  reading. 

3.  Enter  appropriate  data  for; 

a.  AOA 

b.  Dynamic  pressure 
!  c.  Planf orm  area 

:  d.  chord  length 
|  e.  current  temp 
i  f.  delta  P 
!  g.  Tunnel  I  AS 

L  F  or  the  first  pass,  select  "Begin  new 
plot"  with  the  toggle  switch.  The  toggle 
switch  is  located  next  to  the  graphs. 

5.  Take  F orce  calculation  by  clicking  arrow 
in  upper  left  comerof  toolbar. 

6.  Subsequent  runs  must  be  toggled  to 
"Continue  plots". 

7.  Remember  to  change  AOA  each  run. 


For  Information  on  specific 
items  of  the  front  panel; 

1.  Right  click  on  die  item. 

2.  Choose  "data  operations". 

3.  Left  click  "Description". 

To  locate  an  item's  position 
on  the  diagram  from  the  panel; 

1.  Right  click  on  item 

2.  Left  click  “Find  Terminal". 

To  Clear  Graphs: 

1.  Flip  the  Toggle  switch 

2.  Graph  will  clear  wife  next 
data  point  taken. 
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The  first  section  of  the  block  diagram  acquires  the  voltage  signals  (waveforms) 
from  the  strain  gauge  and  determines  the  average  value  for  each  channel.  Output  is  six 
mean  values  for  N 1  through  RM.  The  following  page  converts  the  voltage  signals  into 
forces  (Ibf)  felt  by  each  channel.  It  includes  the  Preinteract  and  Interactl  sub  Vis. 

For  Information  on  a  particular  item: 

1.  Right  click  on  item 

2.  If  necessary  choose  data  operations 

3.  Left  click  on  description  or 

4.  Left  click  on  online  heln 
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Force  Interaction  Equations 
(Interact!  subVI) 


|  Outputs  To 
Tare  And 
Force 

ICalculationsI 


XNl  "EN1 +0.026 7391!2-(2.1944"1  D'-5)Wl2 
ij"l+0.Q0J3233rAX+(1.5099*lQ*-4)'AX"2 
~J-00076122*31+(6i77*10*-3)'S:*2 
-ojocrs?«XKrs2A2 

+0  006«74'RM<37498*10a-5),RM'2; 


XN2'EN2+0.M751*N1-t(66782n0AJ)^r2 
+0002Q7trAX+(33t?ncr^),AX'2 
-O.OOM771-S1-K33223TIO^-5)^S1''2 
-OOC38C36*$2+{1  7359*10A-^),S2A2 
O.00U403"RM-{3i982*10*-4)'RM*2, 


4 


xsi-ESi+ooi3097*m-(umo*j)rNr2 
■0D02123raM(7.771 7KT-9JW2 
+O.OOM761*AX4^  69:H0*-J)*AX*2 
•4Oj055667-S2^-|  1523-10a-4)TS2a2 
■  ■ 000406M*RM-(U997*10*-$)f'RM*2, 


> 


|XAX“EAX+0  0039247*N1+(2.2967’,10*-S)*j'n*2 

1+001 4644"N2+(6.0243*10*'3)?N2''2 

1+0.01 0915*Sl+(6.8636*10*.JJ,sr2 

LoOO*$2+0.00*S2*2 
1+0.01 9642*RM-(1.6028*10* -4  )TRM*2; 


m 

i 

mt-ERMO.OOJ17yNl<4 103*1 0*3)141*2 

43004587*N2-0.0002332*N2*2 

+0.011  74*AX+£!  0762*1 0*-5  FAX*2 

-0.00*S1-0.00"$1*2 

-000*32+000*32*2, 

1 

1 

79 


Results  of  individual  force  calculations  are  routed  to  the  Tare  case  if  the  Tare 
button  is  selected  on  the  front  panel.  Otherwise  they  go  to  the  diagram  on  the  next  page. 
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Remaining  calculations,  including  plotting  and  appending  to  spreadsheet  file,  take 
place  in  this  portion  of  the  block  diagram.  Due  to  the  amount  of  information  displayed,  it 
has  been  split  into  three  sections  in  the  three  pages  that  follow. 


■rrn 
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APPENDIX  B:  STING  BALANCE  CALIBRATION 


A.  Sting  Balance  Calibration  Constants 


Balance  Calibration 

Performed  by:  NASA  Ames  Calibration  Laboratory 


Cal  Date:  8283 

Comp.  Date:  8313 
Invoice  No.:  440528 
Type:  Force 


Size  0.75 

Make  Task/ Able 

pin  No  1 

Rig  No.:  2 


Gage 

Maximum 

load 

Ohms 

X  Gage 

Calibration 

Reading 

N1 

25.00 

25.00 

350 

0.0854 

80K 

5725 

N2 

25.00 

25.00 

350 

0.0854 

80K 

5738 

AX 

50.00 

50.00 

350 

N/A 

80K 

5752 

SI 

12.50 

12.00" 

350 

0.0698 

80K 

5732 

.  S2 

12.50 

12.00 

350 

0.0698 

80K 

5736 

RM 

25.00 

25.00 

350 

N/A 

70K 

6575 

Gage 

KPos(l) 

K  Pos(2) 

K  Neg(l) 

KNeg(2) 

Maximum 

Deviation 

% 

ACC 

N1 

3.6584E-03 

-1.8579E-09 

3.6748E-03 

-1.8922E-09 

-0.054 

•  N2 

3.3823E-03 

-2.6264E-09 

3.3785E-03 

-5.5067E-09 

0.191 

AX 

8.0389E-03 

-7.8379E-09 

8.0496E-03 

0.054 

SI 

1.7347E-03 

-1.6500E-09 

1.7463E-03 

-1.3754E-09 

-0.021 

0.171 

S2 

1.6224E-03 

-1.1579E-09 

—rxmrwri— i 

0.188 

RM 

3.1732E-03 

1.0253E-09 

3.1600E-03 

-1.2994E-09 

0.032 

0.128 

Degree  of  Fit:  2  Accuracy:  15  Int.  Degree  of  Fit:  2 
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Interaction  Coefficients 


Coefficient 

Positive  Value 

N1/N2 

-2.673 9E-02 

Nl/AX 

-5.3233E-03 

Nl/Sl 

7.6  122E-03 

N1/S2 

O.OOOOE+OO 

Nl/RM 

-6.4174E-03 

N2/N1 

-4.3853E-02 

N2/AX 

•  6.2273E-03 

N2/S1 

4.465  8E-03 

N2/S2 

-4.8747E-03 

N2/RM 

-6.1  509E-03 

AX/N1 

O.OOOOE+OO 

AX/N2 

-9.2497E-03 

AX/SI 

O.OOOOE+OO 

•  AX/S2 

O.OOOOE+OO 

A»RM 

-2.1O34E-02 

Sl/Nl 

-1.093 9E-02 

S1/N2 

1.5093E-03 

Sl/AX 

-1.7751E-02 

S1/S2 

-5.5009E-02 

Sl/RM 

9.0376E-03 

S2/N1 

5.5317E-03 

S2/N2 

-2.795 8E-04. 

S2/AX 

8.4602E-03 

S2/S1 

-6.7550E-02 

S2/RM 

6.4729E-03 

RM/N1 

2.2569E-03 

RM/N2 

O.OOOOE+OO 

RM/AX 

-8.4466E-03 

RM/S1 

O.OOOOE+OO 

-2.8717E-02 


-2.2093E-03 


6.9447E-03 


RM/S2 


0.0000E+00 


0.0000E+00 


-4.73  1OE-02 


-2.O761E-03 


3.4771E-03 


3.8036E-03 


4.4405E-03 


-5.9247E-03 


-1.4644E-02 


-1.0915E-02 


O.OOOOE+OO 


-1.9642E-02 


-1.3097E-02 


2.1238E-03 


-3.3761E-03 


-5.5867E-02 


4.0654E-03 


5.18O5E-03 


0.0000E+00 


1  .4573E-03 


-7.2292E-02 


2.761 8E-03 


5.1750E-03 


4.5870E-03 


-1.174OE-02 


O.OOOOE+OO 


0.0000E+00 


Nl/AX 


praam 

i msm 


N2/AX 


N2/S1 


-1.5099E-04 


Coefficient  Positive  Value  Negative  Value 


3.1064E-05 


-8.1729E-05 


2.0633E-05 


0.0000E+00 


O.OOOOE+OO 


-6.6782E-05 


-3.3493E-05 


-3.3223E-05 


-1.7559E-04 


3.6982E-04 


-2.2967E-05 


-6.0254E-05 


0.0000E+00 


3.7498E-05 


1.7830E-05 


1.7338E-04 


-7.499  IE-05 


-5.865  IE-05 


3.7458E-05 


0.0000E+00 


-1.4649E-05 


0.0000E+00 


0.0000E+00 


1.9713E-04 


1.6053E-05 


-1.0662E-05 


4.4896E-04 


8.3835E-07 


-9.2627E-05 


-2.481  IE-05 


-2.0669E-05 


-1.9659E-04 


-1.2240E-04 


-3.3908E-05 


7.9926E-06 


0.0000E+00 


0.0000E+00 


1.6028E-04 


1.6790E-05 


-7.7917E-08 


-5.695  IE-05 


1.1523E-04 


1.1997E-05 


-9.7624E-06 


0.0000E+00 


1.6275E-05 


-3.6922E-05 


9.3172E-07 


4.1030E-05 


2.3320E-04 


-3.0762E-05 


0.0000E+00 


0.0000E+00 


B.  Sting  Balance  Equations 


The  following  equations  perform  the  direct  force  nonlinear  corrections  and  the  force  interaction 
calculations  using  the  coefficients  listed  in  the  previous  section;.  The  first  group  is  used  by  the 
preinteract.vi  while  the  second,  larger  group  consists  of  the  equations  used  by  the  interact 1: vi. 


1 .  DIRECT  NONLINEAR  FORCE  EQUATIONS 

a.  POSITIVE  EQUATIONS 

•  EN 1  =  .0036584  *  N 1  -  1 .8579E-09  *  (N 1  A  2) 
EN2  =  .0033823  *  N2  -  2.6264E-09  *  (N2  A  2) 
EA  =  .0080389  *  A  -  7.8379E-09  *  (A  A  2) 
ESI  =  .0017347  *  SI  -  1.65E-09  *  (SI  A2) 
ES2  =  . 0016224  *  S2  -  1 . 1 579E-09  *  (S2  A  2) 
ER  =  . 003 1732  *  R  +  1.0253E-09  *  (R  A  2) 


b.  NEGATIVE  EQUATIONS 

EN1  =  .0036748  *  N1  -  1.8922E-09  *  (N1  A  2) 
EN2  =  . 0033785  *  N2  -  5.5067E-09  *  (N2  A  2) 
EA  =  .0080496*  A  -  8.8894E-09  *  (A  A  2) 

ESI  =  .0017463  *  SI  -  1 .3754E-09  *  (S 1  A2) 
ES2  =  .0016353  *  S2  -  8.5641E-10  *  (S2  A  2) 
ER  =  .003 1 6  *  R  -  1 .2994E-09  *  (R  A  2) 


2.  FORCE  INTERACTION  EQUATIONS 


a.  POSITIVE  EQUATIONS 

XN1  =  EN1  +  .026739  *  N2  -  2.1944E-05  *  (N2  A  2) 
+  .0053233  *  A+  1.5099E-04*  (A  A2) 

-  .0076122  *  SI  +  6.677E-05  *  (SI  A  2) 

-0  *  S2-0  *  (S2  A 2) 

+  6.4174E-03*  R  -  3.7498E-05  *  (R  A  2) 

XN2  =  EN2  +  .043853  *  N1  -  1.783E-05  *  (NI  A  2) 

-  .0062273  *  A  -  1.7338E-04  *  (A  A  2) 

-  .0044658  *  SI  +7.499 IE-05  *  (SI  A2) 

+  .0048747  *  S2  +  5.865 1 E-05  *  (S2  A  2) 

+  .0061509  *  R  -  3.7458E-05  *  (R  A  2) 

XA  =  EA  +  0  *Nl-0  *(N1a2) 

+  .0092497  *  N2  +  1 .4649E-05  *  (N2  A  2) 
-0  *  SI  -0  *  (SI  A2) 

-  0  *  S2  -  0  *  (S2  A  2) 

+  .021034  *  R  -  1.9713E-04  *  (R  A  2) 
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XS1  =  ESI  +  .010939  *  N1  -  1.6053E-05  *  (N1  A  2) 

-  .0015093  *  N2  +  1.0662E-05  *  (N2  A2) 
+  .017751  *  A  -  4.4896E-04  *  (A  A  2) 

+  .055009  *  S2  -  8.3835E-07  *  (S2  A  2) 

-  .0090376  *  R  +  9.2627E-05  *  (R  A  2) 


XS2  -  ES2  -  .00553 17  *  N 1  +  2.48 1 1 E-05  *  (N 1  A  2) 
+  2.7958E-04  *  N2  +  2.0669E-05  *  (N2  A  2) 
-.0084602  *  A  +  1.9659E-04  *  (A  A  2) 

+  .06755  *  SI  +  .0001224  *  (SI  A2) 

-  .0064729  *  R  +  3.3908E-05  *  (R  A  2) 

XR  =  ER  -  .0022569  *  N1  -  7.9926E-06  *  (N1  A  2) 

-  0  *  N2  -  0  *  (N2  A  2) 

+  .0084466  *  A  +  8.8297E-05  *  (A  A  2) 

-0  *  SI  - 0  *  (SI  A 2) 

-  0  *  S2  -  0  *  (S2  A  2) 


b.  NEGATIVE  EQUATIONS 


XN1  =  EN1  +  .028717  *  N2  -  3.1064E-05  *  (N2  A  2) 
+  .0022093  *  A  +  8.1 729E-05  *  (A  A  2) 

-  .0069447  *  SI  -  2.0633E-05  *  (SI  A  2) 
-0  *  S2-0  *  (S2  A  2) 

-0  *  R-0  *  (RA2)  . 


XN2  =  EN2  +  .0473 1  *  N 1  +  6.6782E-05  *  (N1  A2) 
+  .0020761  *  A  +  3.3493E-05*(AA2) 

-  .0034771  *  SI  +  3.3223E-05  *  (SI  A2) 

•  -  .0038036  *  S2  +  1.7559E-04  *  (S2  A  2) 

-  .0044405  *  R  -  3.6982E-04  *  (R  A  2) 

XA  =  EA  +  .0059247  *  N1  +  2.2967E-05  *  (N1  A  2) 
+  .014644  *  N2  +  6.0245E-05  *  (N2  A  2) 

+  .010915  *  SI  +  6.8636E-05  *  (SI  A2) 

-0  *  S2-0  *  (S2A2) 

+  .019642  *  R  -  1.6028E-04  *  (R  A  2) 

XS1  =  ESI  +  .013097  *  N1  -  1.679E-05  *  (N1  A  2) 

-  .0021238  *  N2  +  7.79I7E-08  *  (N2  A  2) 

+  .0033761  *  A  +  5.6951  E-05  *  (A  A  2) 

+  .055867  *.S2  -  1.1523E-G4  *  (S2  A2) 

-  .0040654  *  R  -  1.1997E-05  *  (R  A  2) 

XS2  -  ES2  -  .0051805  *  N1  +  9.7624E-06  *  (N1  A  2) 

-  0  *  N2  -  0  *  (N2  A  2) 

-.0014573  *  A  -  1.6275E-05  *  (A  A  2) 

+  .072292  *  SI  +  3.6922E-05  *  (SI  A  2) 

-  .0027618  *  R  -  9.3172E-07  *  (R  A  2) 

XR  =  ER-. 005175  *  N1  -  4.103E-05  *  (N1  A2) 

-  .004587  *  N2  -  .0002332  *  (N2  A  2) 

+  .01 174  *  A  +  3.0762E-05  *  (A  A  2) 

-0  *  SI  -0  *  (SI  A2) 

-  0  *  S2  -  0  *  (S2  A  2) 
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c. 


Outputs  using  NASA  Ames  calibration  constants  only 


Table  A.1 


Weight 

load 

point 

N1  POS 

N2  POS 

N 

TOTAL 

N1  NEG 

N2  NEG 

N 

TOTAL 

1.212 

Front 

1.226 

-0.028 

1.198 

-1.32 

-0.101 

-1.421 

1.212 

Middle 

0.59 

0.59 

1.18 

-0.772 

-0.696 

-1.468 

1.212 

Back 

0 

1.194 

1.194 

-  -0.722 

-0.696 

-1.418 

2.212 

Front 

2.22 

-0.011 

2.209 

-2.375 

-0.076 

-2.451 

2.212 

Middle 

1.079 

1.094 

2.173 

-1.25 

-1.18 

-2.43 

2.212 

Back 

-0.023 

2.212 

2.189 

-0.12 

-2.409 

4.217 

Front 

4.219 

0.013 

4.232 

-4.391 

-4.476 

4.217 

Middle 

2.094 

2.05 

4.144 

-2.256 

4.217 

Back 

0.04 

4.155 

4.195 

-0.135 

10.228 

Front 

10.46 

-0.112 

10.348 

-10.239 

-0.279 

-10.518 

10.228 

Middle 

5.1 

4.95 

10.05 

-5.28 

-5.196 

-10.476 

10.228 

Back 

0.018 

10.221 

10.239 

-0.07 

-10.386 

-10.456 

HI 

\wmm 

Weight 

load 

point 

SI  POS 

S2  POS 

s 

TOTAL 

SI  NEG 

S2NEG 

S 

TOTAL 

1.212 

Front 

1.327 

0.08 

1.407 

-1.201 

-0.012 

-1.213 

1.212 

Middle 

0.713 

■HU 

-0.583 

-0.625 

-1.208 

1.212 

Back 

0.106 

1.288 

1.394 

\mmsm 

-1.215 

-1.203 

2.212 

Front 

2.342 

0.074 

2.416 

-2.214 

-0.008 

-2.222 

2.212 

Middle 

1.235 

iilbiECT^ 

2.38 

-1.104 

-1.107 

iilMI 

2.212 

Back 

0.141 

2.384 

0 

-2.198 

BEEI 

4.217 

Front 

4.381 

4.429 

-4.148 

4.217 

Middle 

2.205 

4.363 

-2.026 

-2.204 

-4.23 

4.217 

Back 

0.212 

4.173 

4.385 

0.023 

-4.243 

-4.22 

10.228 

Front 

10.366 

0.128 

10.494 

-10.002 

-0.267 

-10.269 

10.228 

Middle 

5.247 

5.094 

10.341 

-5.071 

-5.203 

-10.274 

10.228 

Back 

0.238 

10.15 

10.388 

0.156 

-10.438 

-10.282 

90 


D.  Local  sting  Balance  Calibration 

Fine-tuning  of  the  Task  Mark  XX  sting  balance  was  difficult  to  achieve.  Previous 
local  calibration  efforts  involved  adjusting  the  individual  excitation  voltages  to  each 
channel.  This  method  was  not  available  as  the  signal  conditioning  board  uses  a  single 
excitation  voltage  for  all  channels.  In  addition,  it  is  not  well  known  how  the  adjusting  the 
voltage  for  one  channel  effects  the  interactions  of  the  other  five  channels.  It  was  possible 
to  mimic  changing  the  voltages  however,  by  adjusting  the  voltage  output  of  each  channel 
by  a  linear  gain  in  Lab  VIEW.  This  could  be  done  before  and/or  after  the  interactions 
were  performed  to  improve  results,  and  each  method  was  attempted.  As  during  the 
verification  of  the  NASA  calibration,  weights  from  one  to  ten  lbs.  were  used. 

The  procedure  chosen  was  somewhat  different  though,  as  forces  were  loaded 
along  the  span  of  the  balance  rather  than  loading  each  individual  channel.  This  was 
accomplished  by  using  the  adapter  sleeve  that  attaches  other  models  to  the  balance. 
Since  it  surrounds  the  balance  with  a  "smooth  friction"  fit,  an  entire  force,  e.g.  Normal, 
Side,  or  Axial,  could  be  loaded  all  at  once.  This  was  thought  to  best  duplicate  how  loads 
would  be  imparted  on  the  balance  by  the  model  during  testing. 

The  balance  was  first  loaded  with  the  adapter  sleeve  and  sling  that  held  the 
weights.  The  weight  of  these  components  was  0.456  lbs.  but  by  using  the  nulling 
potentiometers  on  the  signal  conditioning  board  the  weight  was  nulled  to  zero.  With  the 
gain  set  to  equal  values  on  all  six  channels,  the  balance  was  then  loaded  in  the  S- 
direction  with  weights  of  one  to  ten  lbs.  A  gentle  finger  tapping  on  the  balance  prior  to 
each  measurement  improved  results  by  minimizing  any  possible  hysteresis  from  the 
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previous  load.  Similar  measurements  were  taken  for  the  S+,  N+,  N-  and  AX+  directions. 
The  final  average  value  between  applied  and  measured  load  differed  by  0.059  lbs.  The 
largest  amount  of  error  occurred  at  the  lower  weights  in  the  N+  and  AX  directions. 

Adjusting  the  gain  on  any  of  the  channels  appeared  to  have  only  negative  effects 
on  the  results.  Small  improvements  could  be  gained  on  the  channel  of  interest,  but  would 
adversely  change  the  values  of  the  other  five  channels.  Another  consequence  was  a 
possible  adverse  effect  on  the  channel  being  adjusted  when  loaded  in  the  opposite 
direction.  In  conclusion,  without  an  actual  laboratory  calibration  it  is  doubtful  that 
overall  results  could  be  improved  from  the  initial  settings.  Table  A.2  presents  a  listing  of 
the  applied  loads  and  resultant  forces  read  in  each  direction  measured. 


t 

Table  A.2 


Weight 

S  + 

S  - 

N  + 

N- 

AX  + 

1.001 

1.020 

1.027 

1.101 

1.079 

wassm 

2.001 

1.992 

2.002 

2.113 

2.004 

3.001 

3.005 

3.050 

3.141 

2.941 

3.140 

4.016 

4.037 

4.172 

4.078 

KS&flH 

5.011 

5.011 

4.992 

5.163 

4.969 

6.012 

6.028 

6.031 

6.164 

5.957 

6.172 

7.012 

7.033 

7.015 

7.151 

7.072 

7.063 

8.017 

8.029 

8.010 

8.106 

8.005 

8.171 

9.017 

9.024 

9.028 

9.150 

9.003 

9.020 

10.018 

10.018 

10.006 

10.120 

10.03 

10.180 

92 


APPENDIX  C:  WIND  TUNNEL  WALL  CORRECTIONS 


A.  WALL  CORRECTION  DESCRIPTIONS 

In  an  effort  to  increase  the  accuracy  of  the  data  collected,  tunnel  wall  corrections 
were  applied  to  the  raw  data.  Applying  tunnel  wall  corrections  is  necessary  because  the 
conditions  under  which  a  model  is  tested  in  a  wind  tunnel  are  not  the  same  as  those  in 
free  air.  The  walls  of  a  closed  test  section  produce  changes  in  flow  patterns  that 
influence  the  aerodynamic  results.  When  applying  these  corrections  the  following  areas 
should  be  considered  [Ref.  21]. 

•  Solid  and  Wake  Blocking 

•  Horizontal  Bouyancy 

•  Streamline  Curvature 

•  Downwash  Corrections 

1.  Solid  Blocking: 

The  presence  of  a  model  in  the  test  section  reduces  the  flow  area  and  thereby 
increases  the  velocity  of  the  air  through  the  test  section  around  the  model.  The  solid 
blocking  velocity  increment  at  the  model  is  much  less  than  that  which  would  be 
calculated  from  a  direct  area  reduction  since  the  streamlines  near  the  tunnel  walls  are 
displaced  more  than  those  near  the  model.  Solid  blocking  causes  an  increase  in  local 
dynamic  pressure  which  in  turn  increases  all  measured  forces  and  moments. 
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2.  Wake  Blocking: 

All  models  tested  in  a  wind  tunnel  will  have  a  wake  behind  them  caused  by  the 
detachment  of  the  models'  boundary  layer.  This  wake  has  a  mean  velocity  lower  than 
that  of  the  freestream.  In  order  to  keep  the  mass  flow  constant  from  the  test  section 
entrance  to  exit  the  velocity  outside  the  wake  must  be  higher  than  the  freestream  velocity. 
This  higher  velocity  results  in  a  lower  wake  pressure  and  hence  a  lower  model  base 
pressure  than  would  occur  in  free  air.  The  wake  blocking  effect  increases  with  wake  size 
hence  is  highest  when  large  flow  separations  occur.  Wake  blocking  also  creates  an 
increase  in  dynamic  pressure  on  the  model. 

3.  Horizontal  Bouyancy: 

Almost  all  wind  tunnels  with  closed  test  sections  have  a  variation  in  static 
pressure  along  the  axis  of  the  test  section  (even  with  no  model  present).  A  thickening  of 
the  wall  boundary  layer  causes  a  decrease  in  the  effective  test  section  flow  area.  Test 
section  pressure  normally  decreases  as  the  end  of  the  test  section  is  approached  and  the 
resulting  pressure  variation  over  the  model  produces  an  additional  drag  force.  The  NPS 
tunnel  test  section  is  slightly  divergent  to  help  counteract  this.  In  addition,  the  amount  of 
horizontal  buoyancy  is  usually  insignificant  for  thin  bodies  like  wings,  and  will  therefore 
not  be  considered. 
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4.  Streamline  Curvature: 

The  natural  free  air  streamline  curvature  that  exists  about  any  foil  or  lifting  body 
is  effected  by  the  ceiling  and  floor  of  the  tunnel.  These  tunnel  surfaces  induce  a 
chordwise  variation  in  upwash  and  the  model  acts  as  if  it  had  extra  camber.  Thus,  tests  of 
a  foil  or  lifting  body  in  a  closed  jet  tunnel  will  show  too  much  lift  as  well  as  too  great  a 
moment  about  the  quarter  chord.  Furthermore,  the  effective  angle  of  attack  will  be 
increased. 

5.  Downwash  Corrections: 

The  effect  of  the  sidewalls  of  the  tunnel  is  only  applicable  to  3-D  models  where 
trailing  vortices  exist.  In  this  case  the  tunnel  walls  induce  a  spanwise  variation  in  upwash 
which  decreases  the  normal  downwash.  Thus,  tests  of  a  wing  or  lifting  body  will  have  too 
little  induced  drag  and  the  wing  will  appear  to  have  a  larger  effective  aspect  ratio  than  it 
would  if  tested  in  free  air.  These  effects  are  small  if  the  span  is  much  less  than  the  tunnel 
width  (if  mounted  sideways,  height)  and  in  any  case  will  be  neglected  here. 

B.  CORRECTIONS  FOR  TWO  DIMENSIONAL  TESTING 

All  uncorrected  coefficients  are  defined  in  terms  of  an  uncorrected  dynamic 
pressure  (qu)  which  is  the  q  that  would  exist  if  no  model  were  present  in  the  test  section. 
Neglecting  horizontal  bouyancy,  these  are  the  same  outputs  calculated  previously  using 
the  LabVIEW  program.  The  coefficients  will  first  be  modified  by  the  solid  and  wake 
blocking  corrections  (Esb  and  SWb). 
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1.  Solid  Blocking: 

The  solid  blocking  correction  is  defined  in  terms  of  the  velocity  increment 
produced  (AV)  and  the  uncorrected  test  section  velocity  (Vu). 


_  AV  _  Model _volume) 


V. 


C 


3/2 


where:  Ki  =  0.74  for  a  wing  spanning  the  tunnel  width. 
Ki  =  0.52  for  a  wing  spanning  the  tunnel  height. 
C  =  tunnel  test  section  area. 


The  volume  of  an  airfoil  is  often  approximated  by  the  relation:  0.7(max 
thickness)(chord)(span) [Ref.  21].  Due  to  its  unique  shape  the  waverider  tested  was 
estimated  as  0.5(max  thickness)(planform  area)  =  .0862  ft3.  Wind  tunnel  geometric  area 
=  8.75  ft2,  hence  esb=  0.001732. 

2.  Wake  Blocking: 

For  a  model  which  spans  the  test  section  width  (wall  to  wall)  the  following 
equation  is  normally  used: 


£wb  “ 


AV  cC. 


du 


V. 


2  h 
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For  a  model  which  spans  the  test  section  height  (which  is  the  setup  used): 


A  V  _cCdu 

^  wb  TT 
"  u 

where:  c  -  model  chord 

h  =  test  section  height 
w  =  test  section  width 
Cdu  =  uncorrected  drag  coefficient 

The  coefficient  for  wake  blocking  is  a  function  of  Cau,  hence  a  function  of  AOA.  Therefore,  a 
different  value  is  obtained  for  each  data  point. 


3.  Aerodynamic  Corrections: 

The  total  blockage  correction  (8)  is  found  by  summing  the  solid  and  wake  blockage 

corrections.  Once  the  blockage  factors  have  been  calculated  the  following  expressions  may  be 
used  to  determine  the  corrected  aerodynamic  performance  parameters,  providing  the  chord  to 
tunnel  height  ratio  (c/h)  is  less  than  0.7.  The  subscript  “u”  denotes  uncorrected  values. 

Velocity:  V  =  Vu  (1  +  s) 

Dynamic  Pressure:  q  =  qu  (1  +  2s) 

t 

Reynolds  Number:  Rceff  =  Reu  (1  +s) 

Drag  Coefficient.  du  ~  ^  ^sb  ~  2 £wf  ) 
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AO  A,  Cl  and  Cm  require  correction  for  streamline  curvature,  introducing  a  new  parameter  (c). 


Curvature  Coeff: 

Angle  of  attack: 

Lift  Coefficient: 

Moment  coefficient: 


(7  = 


71 

48 


f-1 

\h) 


(h  =  test  section  height) 


57.3(7/  \ 

a  =  a+ - C,  +Cm 

2tc  '  ci W 

Ci  =  Clu  (\  -  (j-2  s) 


CmM  -  l-2e)  + 


oC, 
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IAS  =  45  mph 


APPENDIX  D:  WIND  TUNNEL  DATA 


del_p 
(cm  H20) 
1.61 


AOA  CL 


IAS  Vel  True 
(mph)  (mph) 
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45  rrph  Bondary  Layer  Corrections 


Corrected  values  =  ’’new/' 


Vnew  37.2 
Renew  4277; 


£^=0.00104 

varies  with  Cd 


CL  new 

CD 

CD  new 

0.0125 

1.2949 

0.7319 

-0.0183 

1.2935 

0.7318 

-0.1118 

1.2703 

0.7284 

-0.1762 

12542 

0.7260 

-0.2394 

1.2177 

0.7196 

-0.3533 

1.2037 

0.7170 

-0.4222 

1.1828 

0.7128 

-0.4660 

1.1647 

0.7089 

-0.5430 

1.0748 

0.6864 

-0.6251 

0.9930 

0.6612 

-0.8216 

1.1317 

0.7013  • 

-0.9748 

1.0423 

0.6769 

-1.0761 

0.8671 

0.6138 

-'1.0889 

0.6467 

0.5053 

-0.9835 

0.4275 

0.3653 

-0.7862 

0.2493 

0.2278 

-0.5460 

0.1214 

0.1161 

-0.2979 

0.0569 

0.0557 

-0.0877 

0.0458 
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0.1608 

0.0715 

0.0696 

0.4431 

0.1524 

0.1442 
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02938 

0.2641 

0.9008 

0.5022 

0.4165 

0.9904 

0.7323 

0.5513 

0.9627 

0.9704 

0.6535 

0.6489 

0.7984 

0.5834 

0.5828 

0.9020 

0.6280 

0.5163 

1.0015 

0.6640 

0.4554 

1.0956 

0.6921 

0.3837 

12235 

0.7207 

0.3239 

1.3219 

0.7353 

0.2579 

1.4093 

0.7429 

0.1938 

1.4734 

0.7452 

0.1319 

0.0697 

1.5619 

0.7438 

0.0105 

1.5571 

0.7441 

-0.0478 

1.5395 

0.7447 

CM 

CM[l/4c) 

-0.2257 

0.4757 

-02434 

0.4588 

-0.2586 

0.4380 

-0.2744 

0.4259 

-02873 

0.4100 

-0.3520 

0.4604 

-0.3598 

0.4375 

-0.3600 

0.4023 

-0.3642 

0.3869 

-0.3699 

0.3805 

-0.5268 

0.4348 

-0.6978 

0.3122 

^•5238228 

-0.6240 

0.2252 

-0.4658 

0.2127 

-0.3080 

0.1867 

-0.1527 

0.1694 

-0.0278 

0.1426 

0.0675 

0.1168 
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60  mph  Boundary  Layer  Corrections 


Corrected  values  =  "new" 
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75  rrph  Tunnel  XAfeill  Corrections  Corrected  values  =  "new/' 
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IAS=  100  rrph 
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lOOrrph  Tinnel  \Afell  Corrections  Corrected  values = "naA/' 
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1 1 0  mph  T unnel  Wall  Corrections  Corrected  values  =  "new" 
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APPENDIX  E:  PRICE  WAVERIDER  PHOTOGRAPHS 


Price  Waverider  As  Seen  From  Outside  Wind  Tunnel 
Turntable  Mechanism  Is  Visible  Below 
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Price  Waverider  As  Seen  From  Settling  Chamber,  Negative  AOA 


Price  Waverider  As  Seen  From  Settling  Chamber,  Positive  AOA 
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